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Abstract 

. Extensions of the Minimal Supersymmetric Standard Model (MSSM) with additional singlet 

(N 
m 
o 

. ing MSSM extensions: Next-to- Minimal Supersymmetric Standard Model (NMSSM), the nearly- 

1^ \ Minimal Supersymmetric Standard Model (nMSSM), and the [/(l)'-extended Minimal Supersym- 
I ' 

Q^' metric Standard Model (UMSSM) by performing scans over model parameters. We find that the 

' Secluded C/(l)'-extended Minimal Supersymmetric Standard Model (sMSSM) is identical to the 



X 



scalar fields solve the important ^u-parameter fine tuning problem of the MSSM. We compute 
and compare the neutral Higgs boson mass spectra, including one-loop corrections, of the follow- 



nMSSM if three of the additional scalars decouple. The dominant part of the one-loop corrections 



H ' are model-independent since the singlet field does not couple to MSSM particles other than the 

Higgs doublets. Thus, model-dependent parameters enter the masses only at tree-level. We apply 
constraints from LEP bounds on the Standard Model and MSSM Higgs boson masses and the 
MSSM chargino mass, the invisible Z decay width, and the Z — Z' mixing angle. Some extended 
models permit a Higgs boson with mass substantially below the SM LEP limit or above theoretical 
limits in the MSSM. Ways to differentiate the models via masses, couplings, decays and production 
of the Higgs bosons are discussed. 
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I. INTRODUCTION 



Supersymmetry (SUSY) is a leading candidate for physics beyond the Standard Model 
(SM). It is the only extension of the bosonic spacetime Poincare symmetry to include a 
fermionic spacetime. Superstring theory, the currently prevailing paradigm of quantum 
gravity, generally includes SUSY, though not necessarily at the weak scale. The cancellation 
of the quadratic divergence in the Higgs mass-squared radiative correction, requiring fine- 
tuning in the SM, strongly motivates SUSY at the TeV scale. TeV-scale SUSY also unites 
the gauge coupling constants at the GUT scale and provides an attractive cold dark matter 
candidate, the lightest neutralino, when R-parity is conserved. The simplest supersymmetric 
extension of the SM is the Minimal Supersymmetric Standard Model (MSSM). 

The MSSM suffers from the yU-problem jl]. The /i-parameter is the only dimensionful 
parameter in the SUSY conserving sector. Naively, in a top down approach, one would 
expect the /i-parameter to be either zero or at the Planck scale, O{10^^) GeV. At tree-level, 
the MSSM gives the relation 



2^ ~ tan2/3-l " ^ ' 



where and m„ are the soft mass parameters for the down-type and up-type Higgs, 
respectively. With the soft parameters at the EW/TeV scale, fi must be at the same scale, 
while LEP constraints on the chargino mass require fi to be non-zero J^. A simple solution 
is to promote the /i-parameter to a dynamical field in extensions of the MSSM that contain 
an additional singlet scalar field that does not interact with MSSM fields other than the 
two Higgs doublets. Extended models thereby circumvent the need for a fine tuning of the 
/i-parameter to the electroweak scale. 

The discovery of Higgs bosons is a primary goal of the Tevatron and the Large Hadron 
Collider (LHC) experiments. Although Higgs boson signals at colliders have been extensively 
studied, most of these studies were based on the assumption that the Higgs bosons occur 
only in doublet fields 141. The few case studies of the Higgs sector in the extensions of the 
MSSM have not been as comprehensive as the SM and MSSM Higgs studies p, 16|. With 
the addition of singlet scalar fields, the properties of the Higgs bosons can be substantially 
different from those in the SM or the MSSM. Moreover, with SUSY, there are also one or 
more extra neutralinos and there may be an extra neutral gauge boson in some models. 
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In this paper we consider models with an extra Higgs singlet field that yield a dynamical 
solution to the /i-problem. The dynamical field that gets a vacuum expectation value (VEV) 
generates an effective /i-parameter that is associated with a new symmetry. These models 
have a third CP-even Higgs boson and, in some cases, an extra CP-odd Higgs boson. The 
mixing with the extra scalar state alters the masses and couplings of the physical Higgs 
bosons. We evaluate the phenomenological consequences of an extra scalar for the Higgs 
masses, couplings, decays and production. We include one-loop radiative corrections to the 
Higgs masses, which to a good approximation turn out to be common among the models at 
this order for the neutral and charged Higgs boson sector. While performing our systematic 
study on the Higgs sector alone, we consider indirect consequences from the neutralino 
sector in anticipation of a later full treatment including both sectors. Detailed studies of the 



neutralino sector in these models have been done by examining the lightest neutralino 
We translate the constraints from LEP experiments on the SM (lightest MSSM) Higgs into 
limits on the CP-even (CP-even and CP-odd) Higgs boson masses in the extended models 
and include constraints from the LEP chargino mass limit, the invisible Z width and the 
Z — Z' mixing angle. 

The extended models of present interest^ are the Next-to-Minimal Supersymmetric Stan- 



oi p 



dard Model (NMSSM) [10|, the Minimal Nonminimal Supersymmetric Standard Model 



(MNSSM) or the nearly Minimal Supersymmetic Standard Model (nMSSM) the t/(l)'- 
extended Minimal Supersymmetric Standard Model (UMSSM) Il2|, and the Secluded t/(l)'- 
extended Minimal Supersymmetric Standard Model (sMSSM) jlS^ . A common /i-generating 
term, hgHu ■ HdS, is contained in the superpotentials of these models, which are listed in 
Table HI After the S field gets a VEV, the effective /i-parameter is identified as 

/ieff = hs{S). (2) 



where {S) denotes the VEV of the singlet field. 

The defining feature of each model is the symmetry that is allowed by the superpotential. 
The NMSSM has a discrete Z3 symmetry, allowing the S'^ term jlfll. With any discrete 
symmetry, the possibility of domain walls exists. It has been shown that domain walls can 
be viewed as a source of dark energy [l^. In the NMSSM, the equation of the state, p = wp, 



^ Many of the ideas of some of the models appeared already in Ref. |3| ■ For a recent review of supersym- 
metric singlet models, see Ref. 



TABLE I: Higgs bosons of the MSSM and several of its extensions. We denote the single CP-odd 
state in the MSSM and UMSSM by A2 for easier comparison with the other models. 



Model 


Symmetry 


Superpotential 


CP-even 


CP-odd 


Charged 


MSSM 




fxHu ■ Hd 






H^ 


NMSSM 


Z3 


hsSHu • Hd + ^S^ 


ttQ ttO ttO 

-"1 ) -"2 ' ^3 


^2 


H^ 


nMSSM 




hsSHu-Hd + CpM^S 


Hi Hi Hi 


4O 4O 
^2 


H^ 


UMSSM 


u{iy 


hgSHu ■ Hd 


Hi Hi HI 


AO 


H^ 


sMSSM 


u{iy 


hgSHu • Hd + XSS1S2S3 


Hi Hi Hi Hi Hi HI 


40 40 40 40 

^2' ^3' ^4 


H^ 



of dark energy is predicted to have w = —2/3 which is disfavored by a recent analysis of 
WMAP data that place w = -1.062+°:J28 

The domain walls may be eliminated if the Z3 symmetry is broken by higher dimensional 
operators, but these may lead to very large destabiliz ing tadpole operators |l7|; one pos- 
sibility for avoiding this problem is described in Ref. 



ig. The nMSSM with a or Zf 



symmetry has a 



domain walls 
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tadpole term of 5* that breaks the discrete symmetries and is thus free from 
|. The harmful tadpole divergences can destabilize the gauge hierarchy, 
but the discrete symmetries or Zy^ allow the divergences to exist only at six and seven- 
loop order, respectively |1S] . At these orders, the divergences are suppressed at scales below 

^Planck- 

An extra U{1) gauge symmetry, U(l)\ is motivated by many models beyond the SM, 



including grand unified theories (GUT) j2Ql, l2l|, ex tra dimensions [2j|, superstrings 
little Higgs [31, dynamical symmetry breaking [2^ and the Stueckelberg mechanism 
The UMSSM and sMSSM each contains a f/(l)' gauge symmetry and its gauge boson, Z', 
that can mix with the SM after symmetries are broken Z f^^. While the continuous 
U(}-)' symmetry is free from domain wall constraints, the UMSSM may require exotic fields 



28l 



29 



30| to cancel chiral anomalies related to the U{iy symmetry^. There are constraints 



on the UMSSM from the strict experimental limits on Z — Z' mixing that are at the mil-level 
js^ . The Z' mass must be above 600-900 GeV to satisfy the Tevatron dilepton search results, 
with the precise experimental limit dependent on the U{iy model With a leptophobic 
Z', these mass limits are evaded. 

^ Exotic fermions can be avoided in a family non-universal U{iy model [s^. 
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The Higgs field content of tlie above listed models is given in Table HJ In the MSSM, 
the usual 2 Higgs doublets give two CP-even -^2)5 ^ CP-odd {A2), and a pair of 

charged (H^) Higgs bosons^. The extended models include additional CP-even Higgs bosons 
and CP-odd Higgs bosons or a Z' gauge boson, depending on the model. The sMSSM 
contains three additional singlets that allow six CP-even and four CP-odd Higgs states. 
However, the additional Higgs fields decouple if A is small and the vacuum expectation 
values (5*1), (5*2), {S3,) are large. The decoupling limit ehminates the D-terms in the mass- 
squared matrix for the S, H^, and fields and yields a model similar to the nMSSM with 
three CP-even and two CP-odd Higgs bosons. This is shown in Appendix ^ We shall 
therefore refer to the nMSSM as n/sMSSM since the results of the nMSSM correspond to 
the sMSSM in the decoupling regime. The charged Higgs sector for all of these models 
remains the same as in the MSSM due to the assumption that the number of Higgs doublets 
is unchanged. 

We present an overview of the Higgs mass-squared matrices including radiative correc- 
tions due to top and stop loops in Section |Tl| We discuss the experimental and theoretical 
constraints applied in Section IIHI and the details of the parameter scans in Section IIVI In 
Sectional we discuss the Higgs spectra and couplings for various models, while implications 
for collider phenomenology are presented in Section IVTl Finally, we summarize our results in 
Section rVlII We provide details of decoupling of the sMSSM in Appendix ^ The derivation 
of the mass-squared matrices of each model are presented in Appendix El and the neutralino 
mass matrices are given in Appendix O In Appendix^ important limits in the Higgs sector 
are addressed, while additional information on the heavier states is given in Appendix lEl 

II. HIGGS MASS MATRICES 
A. Tree-level 

The tree-level Higgs mass-squared matrices are found from the potential, V, which is a 
sum of the F-term, D-teim and soft-terms in the lagrangian, as follows. 

Vf = \h,H^ ■ Ha + ^fMI + kS^\^ + \KS\^ + , (3) 

^ We ignore the possibility of CP-violating mixing effects. 
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Vd = y d^-^l' - + f WlHul' - \Hu ■ H,\') , (4) 

+ ^-^{QHjH,\' + QHjHX + Qsm' (5) 
Koft = ml\H,\' + ml\H^\^ + ml\S\^ + (^A,h,SH^ ■ H, + '^A^S'' + ^sM^S + /i.e.) . 

Here, the two Higgs doublets with hypercharge Y = —1/2 and Y = +1/2, respectively, are 



7- 

















and Hu ■ Hd = eijH^H^. For a particular model, the parameters in V are understood to be 
turned-off appropriately according to Table HI 

NMSSM : gr = 0, = 0, 

nMSSM : gy = 0,k = 0, = 0, (7) 
UMSSM : M„ = 0, K = 0, A, = 0. 

The couplings gi,g2, and gi' are for the U{1)y, SU{2)l, and f/(l)' gauge symmetries, re- 
spectively, and the parameter G is defined as G"^ = gf + The NMSSM model-dependent 
parameters are k and A^ while the free nMSSM parameters are C,f and C,s with Mn being fixed 
near the SUSY scale. The model dependence of the UMSSM is expressed by the D-term that 
has the f/(l)' charges of the Higgs fields, Qh^^ Qhu and Qs- In general, these charges are free 
parameters with the restriction^ that Q ~^ Q Hu ~^ Q s = ^ preserve gauge invariance. In 
any particular f/(l)' construction, the charges have specified values. We assume the charges 
of a Eq model that breaks via the chain Eq — > 5*0(10) x t/(l)^ SU{5) x U{1)^ x f/(l)^ 
[2^. At some high energy scale, the U{1)^ x f/(l)^ symmetry is assumed to break into one 
U{iy^. The above breaking scenario results in the charges 



^ Additional restrictions on the charges of the ordinary and exotic particles come from the cancellation of 
anomalies. 

^ This is the same breaking scheme as in the Exceptional Supersymmetric Standard Model (ESSM) j3J|. 
In the ESSM, among three pairs of SU(2) doublet scalars with MSSM Higgs quantum numbers and three 
singlet scalars, only one pair of doublets and one singlet develop VEVs due to an extra symmetry 
and imposed hierarchical structure of the Yukawa interactions, yielding a model similar to the UMSSM. 
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where OEf^ is the mixing angle between the two f/(l)s and is the only model-dependent 
parameter. 

The F-term and the soft terms contain the model dependence of the NMSSM and 
n/sMSSM. The soft terms of the NMSSM and of the n/sMSSM are new to l/^oft- 

The 5-term of the MSSM is expressed in V^oft as AshgSHu ■ after we identify 



(9) 



The other terms in V^oft are the usual MSSM soft mass terms. 

The minimum of the potential is found explicitly using the minimization conditions found 
in Appendix El The conditions found allow us to express the soft mass parameters in terms 
of the VEVs of the Higgs fields. At the minimum of the potential, the Higgs fields are 
expanded as 

H^a = ^{vd + 4>d + iy^d), H^^ = ^{vu + (pu + iVu), S = ^{s + a + zO- (10) 

with = t>J + t>^ = (246 GeV)^ and tan/5 = Vu/vd- We write the Higgs mass-squared 
matrix in a compact form that includes all the extended models under consideration. The 
CP-even tree-level matrix elements in the S basis are: 



{Ml),, 
[Ml),, 
{Ml),, 
[Ml),, 
{Ml),, 

{Ml) 



2 / hsAs hsKS hs^FMl \ VuS 



(11) 



q2 



[hi + QHaQsgv"^] VdS - 



hsAs hsKs hsiFMl\ 
VdVu - I + H 1 s, (12) 



hgAg 

"Tf" 



V2 



+ hgi^s I V, 



2 , hsAs hsKs hsjpMl \ VdS 
' 2 + s ) v^-' 



[hi + Qh^QsQv'^] VuS - 



hgAg 



+ hsKS Vd 



33 



[n2 2 , o 21 2 , ( hsA, V2^sM^\ VdVu kA^ 
[Qs9i' +2k, \s +1 I \ —s. 



V2 



VdVu 



V2 



(13) 
(14) 
(15) 

(16) 



The tree-level CP-odd matrix elements are: 
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hsAs hsKs h^^FMl\ VuS_ 
V2 2 s ) Vd'' 



/ . .n \ / hgAs hgKs hgfpM^ 

(A^° )^2 = f + + I 



V2 



(17) 
(18) 
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(•^-)23 = (^-^^'^^)^^' (21) 

(A4°)33 = + (22) 

The tree-level charged Higgs mass-squared matrix elements are: 

(-M^)n = ""^^'^^^'^ + [-^^^h,^ + + hsipMl^ ^, (23) 

(-M^),, = _ + \Kns' + /..e^M„^) , (24) 

[M^),, = ""^^^^'^^^'^ + (^^^^^^ + \hs>^^' + ^^e^K') ^- (25) 

The physical Higgs boson masses are found by diagonalizing the mass-squared matrices, 
M.£) = RM-R^^i where M. also includes the radiative corrections discussed below. The 
rotation matrices for the diagonalization of the CP-even and CP-odd mass-squared matrices, 
and for the charged Higgs matrix, T^*-' , may then be used to construct the physical Higgs 
fields. 

Hi = Rl(t)d + B^l^u + R+a, (26) 
Ai = i?!_Vd + R'-Vu + R'X (27) 
= n'^R- + 7^'2i^+. (28) 

where the physical states are ordered by their mass as M^^ < M^^ < Mh.^ and M^^ < ■ 
Many features of the models are apparent by inspection of the mass-squared matrix elements. 
We discuss these aspects in Section El 



B. Radiative Corrections 



An accurate analysis of the Higgs masses requires loop corrections. The dominant con- 
tributions at one-loop are from the top and scalar top loops due to their large Yukawa 
coupling. In the UMSSM, the gauge couplings are small compared to the top quark Yukawa 
coupling so the one-loop gauge contributions can be dropped. Corrections unique to the 
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NMSSM and n/sMSSM begin only at the two-loop level. Thus all contributions that are 
model-dependent do not contribute significantly at one-loop order and the usual one-loop 
SUSY top and stop loops are universal in these models. A similar approach has been done 
in studying extended Higgs sectors with many additional singlet fields js^. These one- loop 
corrections to the potential can be found from the Coleman- Weinberg potential and are 
reviewed in Appendix IB 21 

The mass squared matrix elements become 



M±=Ml+Mi, 



(29) 



where the radiative corrections to the CP-even mass-squared matrix elements are given by 




[m 



2\2 



(m2 — m- Y 



~ 2 ~ 2 

m^m2 



^ , 9 / hgh'jAt ^\ VuS 



hsh^At 
2V2 



T\ s 



~ 2 ~ 2 
■mirrig 



irrY — m~ Y 



hsh'jAt 
2V2 



[m. 



2\2 



m~ — m~ Y 

tl 12' 

hshfAt 



g + 



2hfmn , ^, ,4, 



m~ -|- m~ 

tl t2 



2 2 

m~ m~ , 

tl t2 \ 2 



m 



2V2 



2™2 



mom 



(nY — mi Y fn~ + fn~ 

^ tl t2' tl t2 



[m 



2\2 



(m~ — m~ Y 

^ tl t2' 



g\s' + 



hshfAt^\ VdVu 



2V2 



2V2 



{2-g)]v^s-{^^^J^\v, 



(30) 
(31) 
(32) 
(33) 



(34) 
(35) 



where k 



(47r)2 



and the loop factors are 



g (m| , m 



tl t2 



m~ 



mt 



log 



t2 



mj 



m 



t2 



J^ = los 



2 2 

mf m- 

tl t2 



g{ml,ml). 



(36) 



Here we have defined 



m„ 



htfics (/iefT - AttanP) , 
hfAt {At - fics cot (3) , 
-^h'fues if^es - A tan p) 



(37) 
(38) 
(39) 
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with Q being the DR renormahzation scale and At is the stop trihnear couphng. 
The corrections to the CP-odd mass-squared matrices are given by 

(-Mi). = ^f^^^K,-|), (40) 
where we identify Vi = Vd-, V2 = v^, and ^3 = s. These one-loop corrections agree with those 



of 



19 



The one-loop corrections to the charged Higgs mass are equivalent to those in the MSSM 
and can be significant for large tan (3. The charged Higgs boson in the MSSM has a tree-level 
mass 

(M^°i)2 = M^ + M2, (41) 
and the extended-MSSM charged Higgs boson mass is 



^M^f^M^ + M}-^^^K ^^^^;;'^^\ (42) 

where My = ^^^2^' tree-level mass of the MSSM CP-odd Higgs boson. The case of 

large My (or Ma in the MSSM) yields a large charged Higgs mass and is consistent with 
the MSSM decoupling limit yielding a SM Higgs sector. Radiative corrections in the MSSM 
shift the mass by 

{M^nlf = ^^ + ^Ml,^ (43) 
V2sm 2p 

where after including tan [3 dependent terms, 5M'^± is given by the leading logarithm result 



of the full one-loop calculation 

= 00JL2 „: . J\ a - M'w —h + ]+oMw] log^f^, (44) 



327r2M2, Vsin2/5cos2/3 ^ {3 cos^ {3 ) 3^ 
where A^^^ = 3 is the number of colors and Msusy is the supersymmetric mass scale, taken to 
be 1 TeV. Model-dependent terms come in at tree-level, giving a charged Higgs mass after 
radiative corrections of 



Mf,. - + M# - ^ + + [iyi-^ + ^Mi. 1 . (45) 



III. CONSTRAINTS 



Both theoretical and experimental constraints are important in ensuring that the models 
are realistic. In the following, we list the constraints that we apply in obtaining the allowed 
Higgs mass spectra. 
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To generate the Higgs boson masses, we scan over the relevant parameters of each model. 
Theoretical constraints eliminate large regions of the parameter space. To avoid solutions 
that contain unstable saddle-points of the potential, we require that the mass-squared eigen- 
values are positive-definite, i.e. M^.,M^^,M^± > 0. We also exclude solutions which give 
m~ < 0. 



A. Direct constraints 



The direct constraints are provided by collider data. Currently, LEP gives the best 
experimental bound on the mass of the SM Higgs boson, h, of 114.4 GeV at 95% C.L Q- 
We translate this to limit the mass of the lightest Higgs boson of the extended models by 
using the ZZh coupling limits from LEP, as reproduced in Fig^ 
consider all SM particle decay modes down to = 12 GeV [38' 
relative to the SM coupling is given by the factor 
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of Section IV B| that 
The ZZHi coupling 



^zzm = {gzzHjglfnT = {R^lcosP + R^sinf^f. (46) 

Since the ZZHi coupling in extended models is reduced by doublet-singlet mixing effects, it 
is possible to have Higgs bosons lighter than the SM bound of 114.4 GeV. The reach of the 
ZZHi coupling limit extends only to 12 GeV, below which we do not enforce this constraint. 
However, this low mass region is well constrained by bounds on Ma and in the MSSM 
discussed below. The LEP bound is also applied to H2 and H^ since a heavier Higgs boson 
may violate the bound even if the lightest does not. 

Another channel of relevance from LEP is Z ^ AiHj with Ai bb and Hj bb. 
Current limits place the lightest possible CP-even and odd MSSM masses at Mhj^ = 92.9 
GeV and = 93.4 GeV, respectively and are calculated assuming maximal stop mixing, 



yielding the most conservative limit on the lightest Higgs masses in the MSSM 42|. An 



^ These limits actually assume standard model branching ratios for the Hi, which are dominantly into bb 
and r+r" in the relevant mass range. As discussed in Section FVI Bl for some of the parameter values 
in the extended models the dominant decays are into (invisible) ncutralinos, or into two light CP-odd 
states, and for those points the constraint in Figure [3^^ does not strictly apply. However, there are also 
quite stringent limits on the invisible Hi decay modes [soj . and (weaker) limits on the decays into two 
CP-odd states which subsequently decay into bb or r+T~ 0, 0|. These have not been given for the 
entire kinematic ranges of interest here, so we will simply take the conservative approach of allowing only 
those points satisfying the ZZHi coupling limit in Figure 
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estimation of the corresponding limits in extended models may be obtained by comparing 
the expected production cross section of the extended-MSSM models at the maximum LEP 
energy, ^/s = 209 GeV, to that of the MSSM At this energy the mass limits of the 
CP-even and CP-odd Higgs bosons provide an upper bound of the cross section at 40 fb. 
In practice, we find that the LEP Z AiHj constraint eliminates a significant fraction 
of the points generated with a low CP-odd Higgs mass. In Fig. IHb, we show cos^(« — /?), 
the prefactor of the ZAh coupling where a is the rotation angle required to diagonalize the 
MSSM CP-even Higgs mass-squared matrix, versus CP-odd Higgs mass for the MSSM. A 
strong ZAh coupling results in an enhanced Ah production cross section. In the extended- 
MSSM models, we calculate the cross section for e"^e~ —>■ AiHj where Ai is the lightest 
nonzero CP-even Higgs for that model. If it is above the calculated LEP limit of 40 fb, the 
generated point fails this constraint. Mixing effects which maximize the ZAiHj coupling 
in the MSSM also result in a lower value of M^j, so that the LEP limit implies a lower 
bound on M^j. With the two complementary limits on the neutral Higgs bosons and the 
charged Higgs mass {Mfj± = 78.6 GeV from LEP 4^), the Higgs sector in the MSSM and 
extended-MSSM models are rather well constrained. 



B. Indirect constraints 

While we focus on the Higgs sector in our analysis, indirect constraints from the neutralino 
and chargino sectors also need to be considered. The lightest chargino mass is currently 
limited by LEP to be M^± > 104 GeV at 95% C.L.|2i]. The chargino masses are determined 
by the diagonalization of 

/ M2 y/2Mw sin p 

M^± = 

yV^Mw cos [3 IJ,eG 

The SU{2)l gaugino mass, M2, that enters the chargino sectors does not have a direct effect 
on the Higgs sector, but the lower bound on M^± does constrain possible parameter values. 

Precision electroweak data also provide an upper bound on the new contributions to the 
invisible Z decay width of 1.9 MeV at 95% C.L.^ Contributions to this decay width include 

^ This is based on the constraint on new physics contributions to the invisible Z width, F"^^ = — 2.65± 1.5 
MeV [45! ■ renormalizing the probabihty distribution to require that the true value is positive. Strictly 




12 



Z A,Hj for Ma, + Mh^ < Mz and Z Z*Hi ffH,, for Mh, < Mz- The decay widths 
are given by 

Tz^A.H, = -Mz\'/'(MljMlMljMi){RllV}-RlR^:^)\ (48) 

where X{x,y) = 1 - - - 2{xy - x - y), xh, =_^EhJMz and yu, = MhJMz and 
where the SM Z decay width is Tz^sm = 2.50 GeV \m. Here we assume massless fermions 
in the Z ffHi decay, which is a good approximation at low M/^-. This decay mode 
complements the ZZHi coupling constraint quite well as it is valid below the reach of the 
LEP limit on C,zzh,- Since the masses of H2, and A2 are typically larger than Mz, we 
only consider Z —>■ HiAi and Z ffHi decay modes®. 

The neutralino sector also provides constraints on the allowed parameter space via Z 
boson decay. If M^o < Mz/2, then Z decays into neutralino pairs and this decay contributes 
to the invisible Z-decay width. Since the Z does not couple to the singlino, the superpartner 
of the Higgs singlet, the decay width formula in the extended models is similar to that of 

n 

the MSSM, except for mixing effects |46||. The Z decay width to neutralino pairs, when 
kinematically accessible, is 

^^-^^ = 1^(1^-1' - 1^-"'^' (^^ - (2M,o)f' ' • (50) 
The neutralino rotation matrix elements, Nij, are found by diagonalizing the model- 
dependent neutralino mass-squared matrices in Appendix O 
The Z — Z' mixing angle, 

1 1 / 2M2 , \ 

is also constrained by electroweak precision data to be less than (9(10^'^), where the exact 
value is dependent on the f/(l)' model. The Z' mass parameters are 



Ml, = 9v\Qlvl + Qiy^ + Qls'), Ml,, = gyJgl + gl {v^Qh, - vIQh.) . (52) 



speaking, such decays may not be invisible, and slightly weaker constraints would be obtained using the 
total or hadronic widths. We use the invisible width to be conservative and for simplicity, since it is also 
applied to decays of the Z into neutralino pairs. 

Singlet mixing may allow H2 or A2 to be slightly less than Mz but the decay is still kinematically 
inaccessible. 
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Eq. bounds what types of Z' models and associated Higgs sector parameters are 

allowed; it translates into a high value of s, typically at the TeV scale. There do, however, 
exist isolated points that allow a suppression of azz' at low s such as the following 



(i) If Qh^iQh^ have the same sign, a cancellation occurs at tan/5 = 

(ii) If Qh^ is small and tan is large, the mixing term is suppressed. 



The Z' mass is also constrained j32l . l33| , but the limits are very mo del- dependent on the 
quark and lepton couplings and can be eliminated entirely in the leptophobic Z' case |47|. 
In any case, the large s limit yields a Z' with mass typically large enough to avoid existing 
experimental constraints. Therefore, we only apply the Z — Z' mixing constraint in our 
study. 

Constraints due to the possibility of electroweak baryogenesis have been previously ex- 
plored in the NMSSM the sMSSM Q and the nMSSM 3- The cubic {A,) term in 
the tree-level potential makes it much easier to achieve the needed strongly first-order phase 
transition in these models than in the MSSM 0|. However, we do not consider CP- violating 
phases in the Higgs sector, which is also a necessary condition for baryogenesis. Further- 
more, there are other possibilities for baryogenesis. Therefore, electroweak baryogenesis 
constraints are not included here. 



IV. PARAMETER SCANS 

To generate the Higgs boson masses, we perform both grid and random scans over the 
allowed available parameter space of each model. In the random scan, we evaluate 500000 
points in the available parameter space for each model. Our grid scan gives a reproducible 
catalogue of the Higgs masses of each model. However, due to the large number of parame- 
ters, a finely spaced grid on individual parameters is not feasible. The results from the grid 
scan serve as a useful guide of the allowed Higgs boson masses but do not provide definitive 
upper or lower mass limits. 

The model-independent parameters scanned over are tan s, /Xefr, ^s, and M2, where 

5„2 

we always assume gaugino mass unification Mi = My = g|^2. The masses Mfj and Mq are 
the soft masses of the up-type squarks and doublet-type squarks, respectively, and are fixed 
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TABLE II: Parameter ranges in scans, (a) Model-independent parameters, (b) Model-dependent 
parameters. Parameters not scanned assume the values = Mfj = 1 TeV and Q = 300 GeV. 





tan P 


s 


Mcff 


As 


At 


M2 


Range 


1, 50 


50 , 2000 GeV 


50 , 1000 GeV 


, 1 TeV 


-1 , 1 TeV 


-500, 500 GeV 


Step size 




100 GeV 


100 GeV 


100 GeV 


250 GeV 


100 GeV 



(a) 







Ak 


^5 


Cf 


0Ee 


Range 


-0.75 , 0.75 


-1 , 1 TeV 


-1 , 1 


-1 , 1 


0,7r 


Step size 


0.25 


250 GeV 


0.2 


0.2 


10 



(b) 



at 1 TeV; Mi/ is the mass of the Z'-ino in the UMSSM. The model-dependent parameters 
are k and for the NMSSM, and is for the n/sMSSM, and OEf, for the UMSSM. 

In the parameter scans, we veto points that fail the direct and indirect constraints of 
Section ITTTl We choose the phase convention As > 0, fics > 0, with all the VEVs real and 
positive. We limit hs to be real and positive and allow the gaugino mass M2 and coupling 
n to be real with either sign, although more generally these parameters could be complex. 
With complex parameters, CP violation could occur. If phases were included, the Higgs 
sector would be further complicated with up to five states for the NMSSM and n/sMSSM 
(four states for the UMSSM) that can intermix. The Higgs sector with an arbitrary number 
of additional singlets and CP violation was studied in Ref. js^. 

The couplings run as the energy scale is varied. Naturalness and the requirement that the 
couplings remain perturbative at the GUT scale limits OA < hg < 0.75 or 0.1 < k"^ + h"^ < 
0.75 for the NMSSM. The couplings in the n/sMSSM are real and fixed in the interval 
— 1 < isj^F < 1 with Mn = 500 GeV. We also constrain /igfr to the range 50 < /igfr < 1000 
GeV to avoid fine tuning. A summary of the scan ranges over model parameters are given 
in Table |n] For the grid scan, the step size for each parameter is given and we specifically 
scan tan /3 = 1, 1.5, 2, 10, 50. 
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V. DISCUSSION OF THE HIGGS MASS SPECTRA 

Throughout most of the parameter space, model distinguishing features are apparent in 
the Higgs masses. However, different models can produce similar masses and mixings in 
certain limits. Characteristics that are a direct consequence of how the singlet states mix 
affect the limits placed on the lightest Higgs boson mass. 

A. Common Characteristics 

If the model-dependent parameters in the Higgs mass-squared matrices are set to zero, 
we obtain common mass-squared matrices and an additional symmetry that applies for each 
model. For the n/sMSSM, this is a Peccei-Quinn (PQ) symmetry which protects the mass 
of one CP-odd Higgs. Depending on what parameters vanish, the NMSSM may either have 
a PQ or a f/(l)i? symmetry js^, the global invariance of supersymmetry. Near these limits, 
the Ai mass in these extended models is small, allowing decay modes involving light CP-odd 
Higgs bosons; this is addressed in more detail in Section IVl Bl 

In the UMSSM in the gir limit, the gauged U{iy turns into a global U{1)pq symmetry 
for the matter fields. A massless CP-odd state, Ai, emerges, which is just the Goldstone 
boson of the broken U{1) while the other CP-odd state, present for gir ^ 0, remains massive. 
The Z' decouples and remains massless in this limit. In Table UTTl we summarize the common 
limits of the extended models. 

In the PQ limits (and for the UMSSM for all gii), the CP-odd Higgs mass-squared matrix 
factors into a tree-level matrix times the one-loop correction. Such a form is required by the 
f/(l) symmetries to require the existence of two massless CP-odd goldstones, one of which 
is eaten by the Z and the second by the Z' in the UMSSM after radiative corrections are 
included. Thus, M4 is elevated by a factor of 1 + where the JF term is the loop 

contribution, i.e.. 




(53) 



Effectively the soft mass is increased by 



A,-^As + k^AtT 



(54) 
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TABLE III: Common Higgs mass-squared matrix limits of various models and their effects. Note 
that in the UMSSM, the U{1) is a global symmetry and not a remaining U{iy symmetry. In 
these limiting cases, two of the CP-even Higgs bosons of each model are equivalent to the MSSM 
Higgs bosons if s ^eflfi while the third decouples and is heavy for the NMSSM, or light for the 
n/sMSSM or UMSSM. 



Model 


Limits 


Symmetry 


Effects 


MSSM 




U{1)PQ 


Ma, ^ 


NMSSM 




U{1)PQ 


Ma, ^0 


NMSSM 




U{1)r 


Ma, ^0 


n/sMSSM 




U{1)PQ 


Ma, ^0 


UMSSM 


51' ^0 


U{1)PQ 


Mz',Ma, ^0 



to promote the tree-level mass of the CP-odd Higgs boson to the radiatively corrected one. 
In the U{1)r limit of the NMSSM, the radiative correction to the CP-odd masses vanishes. 

Another limit, the s-decoupling limit, s oo while keeping n^s = ^ ~ C(EW), gives 
similar EW/TeV scale Higgs boson masses for all models. In this limit there is little mixing 
among Higgs states. For the NMSSM and UMSSM, two CP-even Higgs correspond to the 
MSSM Higgs states, while the remaining Higgs boson is dominantly singlet with the mass 
ordering depending on As,k,Ai^ and gr. In the n/sMSSM, the lightest Higgs boson has 
vanishing mass and is singlet dominated while H2 and correspond to the MSSM Higgs 
bosons. Mass expressions in this limit cay be found in Appendix O The Higgs boson that 
is dominantly singlet couples weakly to MSSM particles. 

The strength of a particular Higgs boson. Hi, coupling to fields in the MSSM may be 
quantified as the MSSM fraction 

2 

E(^+)'- (55) 

i=i 

This quantity is not to be confused with the scaled ZZHi coupling ^zzHi- Since R is unitary, 
a sum rule exists 

3 

^MSSM = 2, (56) 

i=l 

which implies that at most two CP-even Higgs bosons can be MSSM-like; equal mixing 

17 



scenarios have ^^^-^'^ ~ §• ^ similar quantity can be found for the CP-odd Higgs bosons. 
In the UMSSM and in the hmits in Table |Tnl for the NMSSM and n/sMSSM, the MSSM 
fraction of the massive CP-odd Higgs boson is 

CmIsm = (l + 5 cos^ sin2 , (57) 
consistent with the s-decoupling limit. 

Since the trace is invariant under rotations, a mass-squared sum rule exists. The limits 
in Table UTTl lead to a common sum rule of the tree- level Higgs masses: 

Tr [Ml - ] = Mlo + Mlo + Ml, - M% = Ml, (58) 

where the Z mass is given by M| = ^(fj + f^). The sum rule for the MSSM is realized 
by taking the s-decoupling limit in the n/sMSSM, and additionally requires gy — in the 
UMSSM and k — in the NMSSM. In the CP-even and CP-odd mass-squared matrices of 
Section m we see that the upper left 2x2 submatrix is that of the MSSM while the third 
column/row vanishes. Then, the decoupled M^^ and M^^ {Mz' for the UMSSM) become 
massless at tree-level and the Higgs mass-squared sum rules become MSSM like: 

M^o + M^o - Ml = M|, (59) 

where hf^ = H2 and = are the usual MSSM CP-even Higgses. 

B. Distinguishing Characteristics 

The introduction of the singlet Higgs field in MSSM extensions produces Higgs boson 
properties that are distinct from those of the MSSM. Each model has additional defining 
characteristics that may be used to distinguish one model from another. In this section, 
we give bounds on the lightest CP-even Higgs mass and provide expressions for the masses 
utilizing the hierarchy of matrix elements given in Appendix |D1 We scan over relevant model 
parameters to determine their effects on the Higgs masses. Finally, we summarize the results 
of the complete random and grid scans. 

1. Lightest CP- even Higgs Mass Bounds 

In any supersymmetric theory that is perturbative at the GUT scale, the lightest Higgs 
boson mass has an upper limit ||52|]. Since the mass-squared CP-even matrix is real and 
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symmetric, an estimation of the upper bound on the smallest mass-squared eigenvalue may 
be obtained by the Rayleigh Quotient 



where u is an arbitrary nonzero vector. With the choices 

= (cos/3, sin/3) [MSSM] 

= (cos P, sin /3, 0) [extended models] (61) 

the well-known upper bound of the lightest Higgs mass-squared from the mass-squared 
matrices of E g. (111111^ and ()30II35|) are given as 
(i) MSSM [5^: 

Mlo < M| cos^ 2(3 + M , (62) 



where 

A^fW = {M^_l%cos^ (3 + {M^_l%sm^ (3 + iM^_^%sm2(3. (63) 

n 

(ii) NMSSM, n/sMSSM, and Peccei-Quinn limits [54]: 

M^o < M| cos^ 2(3 + sin^ 2(3 + M^^\ (64) 

I — I 

(iii) UMSSM [55]: 

M^o < M| cos' 2(3 + hiy sin' 2(3 + gWiQn, cos' /3 + Qh^ sin' /3)' + M^^\ (65) 

Although the upper bounds change with the choice of the u vector, these results indicate 
that extended models have larger upper bounds for the lightest Higgs due to the contribution 
of the singlet scalar. The UMSSM can have the largest upper bound due to the quartic 
coupling contribution from the additional gauge coupling term, gii, in the f/(l)' extension. 
In the MSSM, large tan (3 values are suggested by the conflict between the experimental 
lower bound and the theoretical upper bound on Mh^- Since the extended models contain 
additional terms which relax the theoretical bound, they allow smaller values for tan (3 than 
the MSSM. 



2. Numerical Evaluation of masses 
a. CP-even Higgs Masses 
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In Fig. n we show the variation of the hghtest Higgs mass in the different models as 
functions of s and tan/3 with the other parameters fixed. (Similar plots for the heavier 
states are shown in Appendix |E|) We only apply the theoretical constraints to these spectra 
to see the general trends of the models before experimental constraints are applied. The 
UMSSM would fail to pass the azz' constraint in most of the plotted range of s. 

Note that the MSSM does not conform to the behavior of the extended models in the 
CP-even sector. Since the MSSM contains only two CP-even Higgs bosons, the heavier of 
the two mass-squares increases with fies^s at tree-level, similar to the CP-odd and charged 
Higgs masses. Since we fix hg = 0.5, this Higgs mass-squared scales as the singlet VEV, s. 
The radiative corrections do not contribute a significant s dependence to the mass-squared 
matrix. The tree-level dependence on s prevents a level crossing between the Hi and H2 
states. However, in the extended models there are three CP-even Higgs bosons. Level 
crossings are possible here as there is a Higgs boson of intermediate mass: see Fig. ^c). We 
also see a significant difference between the MSSM and the extended-MSSM models in the 
tan (3 scan, which is expected since a moderate value of s = 500 GeV is chosen. The terms 
that differentiate the matrix elements in the extended models from that of the MSSM are 
not negligible at this value of s, giving different s-dependences of the Higgs mass. 

The MSSM tan (3 scan shows a dip in the Higgs mass at tan (3=1 and a maximal mass 
is approached as tan/5 increases. However, the extended-MSSM models have a decrease 
in mass after tan /3 of 2 — 4 due to the level crossing with the additional moderate mass 
CP-even Higgs present in these models. The presence of the dip in the masses at tan /5 ~ 1 
for the UMSSM and n/sMSSM is not a consequence of a level crossing, but is due to the 
mass dependence on tan (3. When tan (3 = 1, the upper bound on the lightest CP-even Higgs 
mass decreases as seen in Eq. (|62II65|) . Overall, we see substantial differences in the spectra 
of the lightest higgs in the extended models compared to the MSSM. 

b. CP- odd Higgs Masses 

Since only one massive CP-odd Higgs boson exists in the MSSM, UMSSM and the Peccei- 
Quinn limit of the extended models, the CP-odd masses generally behave the same over both 
scans and conform to the general scaling ~ 74ss(cot (3 + tan/5). (The exact expression 
in these cases is given by Eq. fl53|) . with the first term omitted for the MSSM.) Further, we 
note that the CP-odd mass in the Peccei-Quinn limit is identical to that of the UMSSM, 
which may be understood by the absence of mixings and the resulting mass sphttings that 
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CP Even CP Odd 




s (GeV) s (GeV) 

(c) (d) 

FIG. 1: Lightest CP-even and lightest CP-odd Higgs masses vs. tan/3 and s for the MSSM, 
NMSSM, n/sMSSM, UMSSM, and the PQ Hmits. Only the theoretical constraints are applied 
with s = 500 GeV (for tan /?- varying curves), tan/3 = 2 (for s- varying curves). Input parameters 
of As = 500 GeV, At = 1 TeV, Mq = M^y = I TeV, k = 0.5, A^ = -250 GeV, = 500 GeV, 
= —0.1, = —0.1, hs = 0.5, 9e& = — tan~^ ^J^, and Q = 300 GeV, the renormalization scale, 
are taken. The U{\)pq limit allows one massive CP-odd Higgs whose mass is equivalent to that 
of the UMSSM CP-odd Higgs. 



occur in the MSSM or other extended models. However, the MSSM mass approaches the 
PQ/UMSSM mass as s increases, a result consistent with the s-decoupling limit. The lightest 
CP-odd Higgs in the n/sMSSM and the NMSSM, however, does not share the similarities 



21 



of the other models. In these models, there are two CP-odd Higgs bosons, resulting in a 
different dependences on s and tan/3. Mixing effects tend to lower the lightest Higgs masses 
in these models, providing interesting phenomenological consequences. These are further 
discussed in Section IVTl 
c. Higgs Mass Ranges 

We summarize the available ranges found in the grid and random scans of the lightest 
CP-even, CP-odd and charged Higgs boson masses that satisfy the applied constraints in 
Fig. 121 For each model, the values of the maximum and minimum masses are given as well 
as the reason for the bounds. 

The lightest CP-even and CP-odd and the charged Higgs boson mass ranges differ signif- 
icantly among the models. The CP-even Higgs mass range is quite restricted in the MSSM 
and satisfies the upper theoretical mass bound and lower experimental bound from LEP 
discussed in Section IIIII The upper limits for the CP-even Higgs masses in the extended 
models saturate the theoretical bounds and are extended by 30 — 40 GeV compared to the 
MSSM while the upper limits in the lightest CP-odd Higgs masses are artificial in the MSSM 
and UMSSM as they change with the size of the scan parameters such as Ag and tan (3. The 
lower limits of the lightest CP-odd masses in the MSSM and UMSSM reflect the LEP limits 
on Maz', the UMSSM is similar to the MSSM since s is required to be large by the strict 
azz' constraint, decoupling the singlet state and recovering a largely MSSM Higgs sector. 
However, fine tuning the Higgs doublet charges under the f/(l)' gauge symmetry and tan/5 
allows the Z — Z' mixing constraint on s to be less severe, and can result in a lower Higgs 
mass with respect to the MSSM. These instances along with the values Ag = At = GeV 
allow very low CP-even Higgs masses at 0{1 GeV) and a massless CP-odd state. Since these 
points are distinct from the range of masses typically found in the UMSSM, we do not show 
these points in Fig. |21but simply note that they exist. However, the NMSSM and n/sMSSM 
may have a massless CP-odd state due to global U{1) symmetries discussed in Section IVAI 
while the upper limit on the lightest CP-odd Higgs mass depends on the specifics of the 
state crossing with the heavier state, A2, that has a scan-dependent mass. In these models, 
the CP-odd masses extend to zero since the mixing of two CP-odd states allow one CP-odd 
Higgs to be completely singlet and avoid the constraints discussed above. 

The charged Higgs masses are found to be as low as 79 GeV in the scans, in agreement 
with the imposed experimental limit of 78.6 GeV. In these cases where Mj{± ~ 80 GeV, 



22 



CP-Even Higgs Mass Range 



CP-Odd Higgs Mass Range 



2 

Scan 

r_""_ 

2 



LEP 
92 



50 



164 

Th. 

riTj 

170 
Th. 
=1 



100 

Higgs Mass (GeV) 



150 



200 



LEP 

f ■ ■ - 


1 


Scan 


94 




1 

6850 


Th. State Crossing 






1632 






Th. State Crossing 
1 






701 






LEP & a^^. 


Scan 




1 


1 




93 


3617 

1 



2000 4000 
Higgs Mass (GeV) 



6000 



Charged Higgs Mass Range 



LEP Scan 

MSSM r ~_~ ~_"7 ~_"7 ~_"7 ~_~ ~_"7 ~_"7 ~_"7 ~ j 

79 6873 
LEP Scan 
83 6841 

LEP Scan 
n/sMSSM ; 

80 4080 

LEP Scan 
UMSSM I I 
80 3590 

ol . 1 . \ . \ 

2000 4000 6000 

Higgs Mass (GeV) 

FIG. 2: Mass ranges of the lightest CP-even and CP-odd and the charged Higgs bosons in each 
extended-MSSM model from the grid and random scans. Explanation of extremal bounds and 
their values are provided for each model. Explanations are Th. - theoretical bound met, value not 
sensitive to limits of the scan parameters; Scan - value sensitive to limits of the scan parameters; 
State Crossing - value has maximum when crossing of states occurs (specifically for Ai and A2 
in the NMSSM and n/sMSSM); LEP - experimental constraints from LEP; azz' - experimental 
constraints in the UMSSM on the Z — Z' mixing angle. 



the charged Higgs is often the lightest member of the Higgs spectrum. However, these cases 
require fine tuning to obtain values of /Xeir > 100 GeV j9|. The upper limit of the charged 
Higgs mass is dependent on the range of the scan parameters as seen in Eq. (|45|) . The 
discrepancy in the upper limit of the charged and CP-odd Higgs mass between the UMSSM 
and MSSM is a consequence of a lower /Xeir in the UMSSM, resulting in a lower My. Large 
values of /lefr are more fine-tuned in the UMSSM than the MSSM since the additional gauge, 



23 



0.01 



0.001 



o MSSM 
a NMSSM 
V n/sMSSM 
UMSSM 
— LEP limit 95% C.L. 



0.4- 



0.2- 



50 100 
M„ (GeV) 



150 



60 80 100 120 140 
(GeV) 




o oo 
160 " 180 "200 



(a) 



(b) 



FIG. 3: (a) LEP limit Q on ^zzH, = {gzzHjgf^hf = ^z^ZHjTfizh^ the scaled ZZH^ 
coupling in new physics, vs. the light Higgs mass. The solid black curve is the observed limit with 
a 95% C. L. Points falling below this curve pass the ZZHi constraint, (b) cos^(/3 — a) vs. 
in the MSSM. The hard cutoff shown by the solid green line at = 93.4 GeV is due to the 
constraint on (j{e'^e~ — > AiHi) discussed in Section [ill Al 



gy, and Higgs, hg, couplings often drive Mfj^ < 0. Consequently, CP-odd and charged Higgs 
masses comparable to the higher MSSM limit are not present in the scan. The upper bound 
on the charged Higgs mass in the NMSSM is relaxed due to the additional parameter of the 
model. 

d. Higgs Boson Searches 

The focus of Higgs searches is most commonly the lightest CP-even Higgs boson. In 
the models that we consider, the lightest CP-even Higgs boson can have different couplings 
than in the SM. In Fig. IHK, we show the present limits from LEP on the scaled ZZHi 
coupling.^ Mixing effects can lower the ZZHi coupling and, in the MSSM, this occurs if 
is low, as seen in Fig. [Hb where the ZZHi coupling is lowest for cos^(/3 — a) = 1. 
However, an additional limit is placed on the mixing via the e^e~ AiHi cross section 
discussed in Section lilIA[ eliminating low mass CP-even Higgs bosons in the MSSM, as seen 
in Fig. Eb- In extended-MSSM models, additional mixing may occur with the singlet fields. 
Due to this mixing and the subsequent evasion of the LEP limit on the ZZHi coupling, 

^ For clarity, in all the plots that follow we sample the passed points in the results from the random scans. 
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FIG. 4: Higgs masses vs. ^mssm in the (a) NMSSM, (b) n/sMSSM, (c) UMSSM and (d) the 
hghtest CP-even Higgs of all extended models. The vertical line is the LEP lower bound on the 
MSSM (SM-hke) Higgs mass. 



the lightest CP-even Higgs may then have a mass smaller than the SM Higgs mass limit. 
Indeed, attempts to explain the 2.3a and 1.7a excess of Higgs events at LEP for masses of 
98 GeV and 114.4 GeV, respectively, with light CP-even Higgs bosons in the UMSSM have 

been explored 5^. This slight excess has also been studied in the NMSSM where a light 

f I 

Higgs with a SM coupling to ZZ decays to CP-odd pairs 57]. 

The reduction in the CP-even Higgs mass in extended models can be seen in Fig. 01 where 
we plot the MSSM fraction versus the Higgs boson mass. When there is little mixing between 
the singlet and doublet Higgs fields, the MSSM limit is reached and the LEP bound applies, 
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FIG. 5: Higgs mass dependence on ^5 in the n/sMSSM. When ^5 ~ —0.1, H2 and Hi switch 
content, allowing a light CP-even Higgs below the LEP limit. 

as seen by the MSSM cutoffs at ^mssm = 1 and Mh^ = 114.4 GeV. A common feature of each 
model is a CP-even Higgs boson with a mass range concentrated just above the LEP SM 
mass hmit shown by the dark-green vertical line. These Higgs bosons have a large MSSM 
fraction, for which the ZZHi coupling limit is effective in elimination of the generated points. 
We note that there are cases where a Higgs boson mass below 114.4 GeV but with relatively 
high MSSM fraction is allowed due to cancellation between the rotation matrices in Eq. 

This cancellation permits the lightest MSSM Higgs boson to be below the SM limit, 
and has been taken as a possible explanation of the Higgs signal excess Q]. 

By measuring the lightest Higgs boson couplings to MSSM fields, an estimation of the 
MSSM fraction may be obtained, providing important information on the singlet content. 
In the NMSSM and especially the n/sMSSM the lightest CP-even Higgs boson may have 
both low MSSM fraction and low mass as seen in Fig. |31i. Since yUefr is fixed at the EW 
scale, the matrix elements {Ai+)i3 are suppressed in the n/sMSSM at large s. This results 
in a low mass CP-even Higgs boson with high singlet composition; the other Higgs states 
have a high MSSM fraction due to the sum rule in Eq. ()56p. However in the n/sMSSM, the 
existence of a low mass CP-even Higgs boson depends on the value of ^5. In appendix |D 2 g 
we show that the tree-level mass-squares of the singlet dominated CP-even and odd Higgs 
bosons in the n/sMSSM at large s are 
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FIG. 6: Higgs masses vs. s in the (a) NMSSM, (b) n/sMSSM, (c) UMSSM and (d) the Hghtest 
CP-even Higgs of all extended models. The vertical line is the LEP lower bound on the mass of 
the SM Higgs. 



which forces the parameter to be negative in this Umit. Therefore, a largely singlet 
CP-even Higgs boson can have a mass lower than the LEP limit if 

(114 GeVfs ^ ^ 

-^s< 7= 0.1. (67) 

V2MI ^ ' 

In Fig. |S[ we show the Higgs mass dependence on this parameter, which exhibits the crossing 
of states at = —0.1. 

In the UMSSM, the lightest Higgs mass is concentrated near the LEP limit with ^mssm 
near one, which is a direct consequence of the high s constraint placed by the strict azz' 
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limit. This is also seen in Fig El where we plot the Higgs masses versus the singlet VEV. The 
lowest allowed point in the UMSSM has s above ~ 800 GeV, compared to the other models 
which allow s to be as low as a few hundred GeV. By examining Fig. ^ and IHfc we see that 
varies linearly with s and is characteristically dominantly singlet. Without the azz' 
constraint, the Hi and H2 states cross near s ~ 400 GeV. This constraint may be evaded by 
the fine tuning cases discussed in Section IV B 21 At this point, the mass eigenstates switch 
content, below which the lightest Higgs is dominantly singlet, has a mass below the LEP 
bound, and evades the ZZHi coupling constraint. 

The Higgs mass dependence on tan/5 has some interesting features, especially that of 
the lightest Higgs. We show this dependence in Fig. [3 for all the Higgs bosons of each 
extended MSSM model and separately for the lightest Higgs in all the models considered. 
The lightest CP-even Higgs boson mass vs. tan/? in each model shown in Fig. [7|i has a 
majority of generated points in the band 114.4 GeV < Mh-^ < 135 GeV and tan/5 > 2. This 
is one of the salient features of the MSSM as shown in Fig. The MSSM parameter space 
has a lower cutoff at tan/? ~ 2 due to the LEP limit at 114.4 GeV for a SM-like Higgs and 
is shown in FigOi as the intersection of the theoretical MSSM Higgs mass limit shown in 
blue and the LEP limit in green. However, the extended-MSSM models may have values of 
tan (3 that are below this region. Since mixing effects can decrease the lightest Higgs mass 
and thereby satisfy the LEP bounds, a strict bound on tan (3 cannot be given. Additionally, 
an increase in the Higgs mass from the MSSM theoretical limit shown in Section FV B II can 
permit low tan/5 scenarios which have masses above the LEP limit. 

Among these models, the heaviest CP-odd Higgs state follows the same dependence on 
tan (3 that was noted above in Section IVB2I The heaviest CP-even Higgs and charged 
Higgs bosons also follow this trend with the charged Higgs boson mass having the same 
tan/? dependence as the CP-odd Higgs mass, see e.g. Eq. fH^ . The heaviest CP-even 
and CP-odd Higgs masses are approximately the same even after radiative corrections. An 
explanation is provided by the mass-squared sum rules that each model obeys, namely 



The sums are over the massive Higgs bosons, and M^mssm is a model-dependent mass 
parameter with values 
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FIG. 7: Higgs masses vs. tan /3 in the (a) NMSSM, (b) n/sMSSM, (c) UMSSM and (d) the hghtest 
CP-even Higgs of all extended models where the blue curve shows the theoretical MSSM mass limit 
with maximal mixing in the stop mass-squared matrix. 
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The term 5M'^ in Eq. ()68|) is due to the radiative corrections, and has a value 

5M^ = Tr [Ml - Ml] , (72) 

that gives an estimate of the effect the radiative corrections have on the Higgs masses. Note 
that the CP-odd radiative corrections, the J-' terms, are cancelled by equivalent terms in 
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the CP-even mass-squared matrix. The radiative corrections alter the sum rule by at most 
(9(100 GeV)^ over most of the scanned range, as seen in Fig. |H1 where we plot the shift 
versus both tan/5 and At. 

The radiative correction contributions to the sum rule are largest for large At and small 
tan /3. Since the top quark Yukawa coupling increases when tan (3 is small, the radiative 
corrections are enhanced at small tan /?, causing larger deviations from the sum rule. Since 
radiative corrections only affect the sum rule by (9(100 GeV), any high mass CP-even Higgs 
boson contribution must be cancelled by a CP-odd Higgs of similar mass. 



VI. COLLIDER PHENOMENOLOGY 



Higgs boson decays are important to consider as they affect signals at colliders. Both 



production and decay modes are relevant in determining whether a given model yie 



ds de- 



tectable Higgs physics. While Higgs searches have been addressed for the NMSSM 
a side-by-side comparison of the NMSSM, n/sMSSM, and UMSSM has not yet been made. 
In the above parameter scans, we calculate the partial decay widths relevant to production 
and branching fraction for decays in these models of various important modes. 
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A. Higgs Production 



62|- We calculate the partial decay 



6l| with the SM Higgs 



At hadron colliders the dominant production of the lightest Higgs boson in the SM 
proceeds through gg fusion and/or Weak Boson Fusion (WBF). The Higgs production cross- 
sections are directly related to Higgs decay widths when the decay channels are kinematically 
accessible at the Higgs mass. Considerable effort has been put into calculating Higgs decays 
beyond leading order 0|. In the SM, MSSM, and NMSSM, numerical codes have been 
implemented to calculate these widths precisely [giI, 
widths of if ^ gg, WW, and ZZ in each model via HDECAY 
couplings modified to reflect the model of interest. Higgs decay to gg occurs via quark 
and squark loops and is calculated at NLO. However, to a good approximation, the loops 
involving heavy squarks are suppressed {s^]- The large squark mass approximation is justified 
for our assumed values of Mq = Mjj = Mj^ = 1 TeV. Therefore, we only consider the SM 
quark loops in the H ^ gg calculation. Decays to weak boson pairs are calculated at 
tree-level as the radiative corrections to the width are negligible. 

In Fig. ini we plot the partial widths of the lightest Higgs boson for the SM, MSSM and 
extended-MSSM models. Since the n/sMSSM and, to a lesser extent, the NMSSM contain 
a very light Higgs with high singlet composition, its decay widths to SM particles are highly 
suppressed. However, from Fig. |3J the second lightest CP-even Higgs boson in the n/sMSSM 
has a high MSSM fraction and often has comparable mass to the lightest Higgs bosons in 
other models. Hence, we also show the decay width of the second lightest Higgs boson in the 
n/sMSSM, as it is characteristically similar to the lightest Higgs boson of the MSSM. The 
decay widths of the lightest Higgs in the MSSM show a large spread with respect to the SM, 
associated with low A2 mass: see Fig. IHb- When ^ Mz, t he masses and couplings of 
the lightest CP-even Higgs approach those of the SM Higgs |64j. In the UMSSM, the azz' 
limit forces the model to be near the s-decoupling limit, resulting in masses, couplings, and 
decay widths that are close to those of the MSSM. Consequently, the UMSSM decay widths 
lie directly on the SM width in Fig. 121 

In the models considered, the Hi mass is typically below the WW and ZZ thresholds. 
Therefore, the off-shell WW* and ZZ* decay widths are evaluated in the MSSM and its 
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FIG. 9: Decay widths for WW*, ZZ*, and gg in the MSSM and extended-MSSM models. Curves 
denote the corresponding SM width. For clarity, not all points generated are shown. 



extensions^" . For the decays of the very hght Higgs boson to occur in the n/sMSSM two 
off-shell gauge bosons are involved, resulting in high kinematic suppression of decay rates. 
In all the models considered, the WW* and ZZ* partial widths are bounded above by those 
of the SM. This is a consequence of the complementarity of the couplings of Hi and H2 to 
gauge fields in the MSSM. The gauge couplings in the MSSM follow the relation 

/ SM \2 _ (MSSM\2 , (MSSM\2 /7q\ 

\yvvh) — {yvvHi ) + {yvvH2 ) ■ y''^) 
More sum rules exist in the MSSM and can be found in Ref. jo^. In extended-MSSM 

In this case the decay width cannot be translated directly into production rates since they require trans- 
verse and longitudinal polarizations of the I^-bosons to be treated separately. However, the gauge coupling 
is equivalent in either case, and its scaling contains the suppression of the production rate. 
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models the gauge couplings are related to the SM couplings by 

9VVH, = 9vvhit^+ cos/3 + sm/3), (74) 

where Qzz^^ = ^^l^^Ow'' ^"^^ 9wwh — W2Mwgnu- The sum rule in Eq. (17^ generalizes to one 
involving three Higgs couplings. Therefore, the coupling of the lightest Higgs boson to weak 
bosons in the MSSM and its extensions is always reduced compared to the SM couphngs. 
LEP constraints require the ZZHi coupling to be below the SM coupling when the Higgs 
mass is below the 114.4 GeV limit. Associated production qq ^ V* VHi for SM Higgs 
bosons can be important at the Tevatron and the LHC for low Higgs masses 0]. The 
corresponding production cross section can be scaled from the SM calculation by the VVHi 
coupling in Eq. (fTIj). 

The Hi gg partial width governs Higgs boson production via gg fusion at hadron 
colliders. We see in Fig. IHl that the gg partial decay width is typically suppressed in 
the n/sMSSM for a low mass of the lightest Higgs boson since it is dominantly singlet. 
However, there is a trade-off in the production cross-section between smaller r{H — > gg) 
and the kinematic enhancement from a lighter Mh^ whose interplay is beyond the scope of 
this paper. The lightest Higgs in the NMSSM and MSSM and the second lightest Higgs 
in the n/sMSSM have decay widths to gg that may be either enhanced or suppressed by a 
few orders of magnitude depending on the Higgs coupling to the internal quarks and their 
interferences. However, the lightest Higgs in the UMSSM and the MSSM in the limit of a 
large CP-odd mass shows no significant deviations from the SM h ^ gg decay width. 



B. Decay Branching Fractions 

Specific decay modes are important for identifying the Higgs boson at colliders. We 
calculate the contributions of bb, cc, ss, t~^t~, WW*, ZZ*, and gg to the total 

decay width of the lightest or MSSM-like Higgs boson in each model using HDECAY after 
modifying the corresponding Higgs couplings. In addition, non-SM decays including XiXi 
and AiAi are calculated since they are often quite light in the extended-MSSM models. The 
decays to 77, and 7Z are also calculated with loops involving quarks, ly^-bosons, charged 
Higgs bosons, charginos, and squarks (which decouple for sufficiently large squark mass). 

The branching fractions of representative decays to the SM particles WW*, ZZ* , bb. 
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FIG. 10: Branching fractions for various modes in the MSSM and extended-MSSM models. Curves 
denote SM branching fractions. 



r+r , 77, and Z7 are presented in Fig. for the lightest CP-even Higgs boson in the 
MSSM, NMSSM and UMSSM, either or H2 in the n/sMSSM, and h in the SM. 

Note that the branching fractions may be larger in the SUSY models than in the SM. 
For instance, in the NMSSM the branching fractions to WW* and ZZ* can be larger than 
the corresponding SM branching fractions, as seen in Fig^J These enhancements are due 
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to the smaller total decay width of the Higgs boson rather than an enhancement of the 
particular partial width and may aid in the H — > W*W* Wjj and Wlu discovery modes 



at the Tevatron 



67|. Since the dominant decay mode is typically to bb in the mass range 



2mh ^ Mh ^ 140 GeV, any decrease in the bb partial width reduces the total Higgs width. 
The Higgs boson couplings to fermions are related to the SM values by 

- 9ffH^, 9uuH. - 9ffk—^, (75) 

where g^^l = Hence, either suppression or enhancement of the partial decay widths 

to fermions is possible, but not to an arbitrary degree, since in the MSSM the rotation 
matrices and tan (3 obey the tree-level relation 

sin2a= t!f ^!^| sin2/3, (76) 

where sin a = R^_l and cos a = i?!^ in the MSSM. A similar expression holds for the extended 
models, which restricts the rotation matrix values. As noted earlier, the coup lings converge 
to SM couplings just as they do for the MSSM in the s-decoupling limit |64|. 

For a SM Higgs boson of mass below 150 GeV, /i — 77 is a significant mode for discovery 
at the LHC The branching fraction for this mode can be enhanced significantly due to the 
modified fermion loops in the n/sMSSM and NMSSM for the same reasons that the H gg 
decay width is enhanced, providing more opportunity for discovery. The Higgs couplings to 
VT-bosons, charginos and also affect the 77 and Z'j branching fractions, shown in Fig. 
ITUl These couplings are reduced from their MSSM values. However, the reduced couplings 
may not necessarily lead to a rate suppression as interference effects can enhance the overall 
partial decay widths. 



1. Non-SM decays 



Decays to non-SM particles can also be important in the extended models. Since the 
lightest neutralino is a dark matter candidate, its production at colliders is of great interest 
to both the particle physics and cosmology communities 6^. We show in Fig. ITTk the 
kinematic region where neutralino production via the decay Hi — XiXi is possible. The 
couplings and masses of the lightest neutralino have been investigated for the models con- 
sidered here 7|, and M^o may be quite small in the n/sMSSM However, in the 
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n/sMSSM most of the kinematic region is disfavored due to a large Xi relic density |7|. This 
is indicated in Fig. ITTk below the red horizontal line at M^o = 30 GeV, which is the lower 
bound of M^o allowed by the dark matter relic density constraint when only the annihila- 
tion through the Z pole is considered. The Z pole is the most relevant channel since the 
Xi in this model is very light (M^o <^ 100 GeV). In principle, other annihilation channels 
such as a very light Higgs may allow the lighter Xi although the pole will be quite narrow 
Goj ]. Furthermore, in the secluded (sMSSM) version of the model, it is possible that the Xi 
considered here actually decays to a still lighter (almost) decoupled neutralino, as discussed 
in Appendix 1X1 
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FIG. 11: (a) Mfji vs. M o in all the models considered. Points falling below the blue line allow 



the decay of the lightest CP-even Higgs to two Xi- (b) Branching fraction of Hi XiXi 



The XiXi partial decay width is given by 

1 



where the Hixlxi coupling is 



iMyMf,^,MyMl] iMl-4M'A \Cj,^^o^o 



+ 
+ 



{giNii - g2Nu + gi'QH^Nie)Nu + V2hsN,3N^, 



il 



(77) 



(78) 



where the expression for the NMSSM in Ref. [3| has been generalized to include the UMSSM 
while the HiXiXi coupling in the n/sMSSM does not contain any mo del- dependence. For a 
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particular model, the irrelevant parameters are understood to be set to zero as in Eq. (jUj). 
The lightest Higgs boson in the n/sMSSM can have a high branching fraction to the lightest 
neutralino as seen in Figs. ITTb . In fact, in this model the XiXi branching fraction can be 
near 100%.^^ This Z decay is seen as missing energy and makes Higgs searches difficult at 
the Tevatron or LHC It has been explored in the MSSM 72] and more generally 

In addition to decays to neutralino pairs, decays involving the lightest CP-odd Higgs 
bosons are relevant in the extended models. In Fig. IT^ we show the possibilities for decays 
involving both Ai and ifj, where Ai is the lightest nonzero CP-odd state for each model. 
The kinematic regions where Z —>■ AiHi and Hj AiAi are given. Even though the Z 
decay is possible in the n/sMSSM and NMSSM, the coupling is suppressed due to the low 
MSSM fraction of both A^ and Hi seen in Fig. EJd. Also shown is the crossing of states 
in the n/sMSSM where H2 and Hi switch content and hence their variation with M^.. 
The lightest CP-even and CP-odd Higgs masses in both the MSSM and n/sMSSM show a 
strong correlation below the LEP limit. In the MSSM, this is evident from Fig. where 
the reduced ZZh coupling occurs when cos^(/3 — a) does not vanish, resulting in a lower 
CP-even Higgs mass. The n/sMSSM correlation is more clearly shown in Fig. El where the 
crossing of states at ^5 ~ —0.1 is discussed. 

The H — > AiAi mode can be significant if allowed kinematically and has been studied 
in the NMSSM {tJ] and in the general singlet extended MSSM via an operator analysis 
Since the lightest Higgs masses are small in the n/sMSSM at low we scan over this 
parameter with a higher density in this region to be near the PQ limit, which gives a lightest 
CP-odd Higgs boson of low mass. In Fig. IT^ . all the points below the line Mh^ = '^Ma^ 
allow this decay; the corresponding partial width is given by 

T{H, - AA,) = Y^^A^^ {mI/M^, Ml/Ml) \Cn,A.A.\\ (79) 

where the HjAiAi coupling, 

Ch.a.a. = Ph.Pa.PaV, (80) 
is determined by the projection operators that parallel the equivalent MSSM operators in 



If we do not assume gaugino mass unification and assume that ficS is fight and My is heavy, then Xi 
is fight and large XiXi branching fractions are possible in the UMSSM, similar to those found in the 
n/sMSSM. For constraints on M^o in the MSSM from supernova data see Ref. [tJ- 
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FIG. 12: (a) vs. M^. showing the kinematics for decays in extended-MSSM models, where 
Ai is the hghtest nonzero CP-odd state for each model. Hi — > AiAi decays are allowed for regions 
below the blue-dashed line. Decays of Z — > HjAi are allowed to the left of the green dark line, 
(b) H AiAi branching fraction vs. Higgs mass. The n/sMSSM parameter ^5 is scanned with a 
higher density at low |^5| to allow low Higgs masses. 



Ref. 



d 



(81) 
(82) 

= 0, where the field 



where we evaluate the potential at the minimum (f)u,d = cr = (pu,d 
values are shifted to the minimum as in Eq. (jlOj) . 

The Hi A2A2 decay is not allowed in the UMSSM. This is because the azz' constraint 
often requires a large value of s resulting in a CP-odd mass above the typical CP-even 
mass, see, e.g., Eq. ()D14|) . However, the AiAi decay is kinematically allowed in both the 
n/sMSSM and NMSSM. When allowed, this decay mode can be dominant in the NMSSM 
as seen in Fig. [T^b- In the n/sMSSM, the Hi AiAi decay mode is suppressed since 
there is no tree-level coupling of three singlet Higgs states in this model and the Hi and 
Ai states are dominantly singlet. However, because Hi is not completely a singlet, these 
decays can still be non-negligible. In Fig. IT^ . we show the Hi AiAi decay width in 
the n/sMSSM versus the product of the MSSM fractions of Hi and Ai. As both MSSM 
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FIG. 13: (a) Hi AxAx decay width in the n/sMSSM vs. the product of the MSSM fractions of 
}i\ and A\. (b) Mai/Mh^ vs. tan /? in the n/sMSSM. The horizontal hne marks the production 
threshold of Hi AiAi. 



fractions vanish, Hi and Ai become singlet dominated, giving a vanishing decay width in 
the n/sMSSM due to the absence of a singlet self-coupling. Nevertheless, this partial decay 
width alone characteristically exceeds the total width of the SM Higgs boson as can be seen 
in Fig. HH 

The MSSM-like second lightest Higgs boson in the n/sMSSM also may have a large 
branching fraction to light A pairs. Since H2 has a large MSSM fraction, the coupling to the 
singlet Ai pairs is not suppressed. In addition, kinematic suppression is absent due to the 
larger mass of H2- In the lightest Higgs boson decay to AiAi in the n/sMSSM, low tan/5 is 
preferred as shown in Fig. IT5b . where the horizontal line marks the production threshold. 
This decay also requires a low Ai mass and results from the near-Peccei-Quinn limit when 
,^5 0. The low tan/3 preference is a result of the larger Hi mass in this region. This 
enhancement is suggested in the one-parameter scans shown in Figs. ^ and El where the 
lightest Higgs mass is peaked at low tan/5 due to the crossing of Higgs states. In contrast, 
while the NMSSM's lightest Higgs mass is maximal at low tan/3, a sharp drop as tan/5 is 
increased is not present, yielding little to no correlation of tan/5 with the existence of this 
decay mode. 

Finally, we show the total decay width of the light CP-even Higgs bosons to SM modes 
and XiXi and A^Ai in the models considered in Fig. HH The total decay width can be 
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FIG. 14: Total decay width for each model. Large enhancements with respect to the SM are largely 
due to the decays to AiAi and XiXi- 

enhanced due to the XiXi ^ind AiAi partial widths. The total width in the MSSM can be 
larger than the SM due to the enhanced couplings of the Higgs to bb when away from the 
MSSM decoupling limit. In the n/sMSSM, Higgs masses above the LEP bound decaying to 
Xi pairs make a contribution to the total width that is no larger than a few MeV. The AiAi 
decays are responsible for the significantly larger total widths in the NMSSM and n/sMSSM. 

VII. CONCLUSIONS 

Extensions of the MSSM that include a singlet scalar field provide a natural solution 
to the undesirable fine-tuning of the /i-parameter needed in the MSSM. After symmetry 
breaking, the singlet Higgs obtains a VEV, generating an effective /i-parameter naturally 
at the EW/TeV scale. While the extensions to the MSSM that we consider each contain 
at least one additional singlet field, S, the symmetries that distinguish each model and 
their resulting superpotential terms provide phenomenologically distinct consequences. We 
made grid and random scans over the parameter space of each model and imposed the LEP 
experimental bounds on the lightest CP-even ZZHi couplings. The limits on M^^ and M^i 
in the MSSM were converted to associated AiHj production cross section limits and imposed. 
We also imposed constraints from the LEP chargino mass limit and new contributions to 
the invisible Z decay width. Within the UMSSM, we enforced an additional constraint on 
the Z' boson mixing with the SM Z. 
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We found the following interesting properties of the considered models: 

(i) The lightest Higgs boson can have a considerable singlet fraction in the n/sMSSM 
and NMSSM. Since the singlet field does not couple to SM fields, the couphngs of the 
hghtest Higgs to MSSM particles are reduced due to the mixing of the singlet field 
with the doublet Higgs bosons, resulting in the e+e" production cross sections being 
significantly smaller. Therefore, in the n/sMSSM and NMSSM, Higgs boson masses 
that are considerably smaller than the LEP bound on the SM Higgs boson mass are 
possible. The upper bound on the lightest CP-even Higgs mass in extended-MSSM 
models is also relaxed due to the contribution of the singlet scalar through the mixing 
of the Higgs doublets and the singlet. The upper limit in parameter scans is increased 
up to 164 GeV for the NMSSM and 170 GeV for the n/sMSSM. The hghtest CP- 
even Higgs mass in the UMSSM can be as large as 173 GeV due to additional gauge 
interactions; however, the lower bound on the lightest CP-even Higgs mass is similar 
to that of the MSSM. 

(ii) A common feature of each model is a CP-even Higgs boson with a mass in a range 
concentrated just above the SM mass limit. At least two CP-even Higgs bosons must 
have nonzero MSSM fractions, making the lightest non-singlet dominated Higgs obey 
limits on the ZZHi coupling, forbidding masses below 114.4 GeV unless additional 
doublet mixing occurs. The lightest Higgs in the NMSSM and n/sMSSM can evade 
this by singlet-doublet mixing, allowing a Higgs with mass just below the SM limit. 

(iii) In the s-decoupling limit with fixed //eS at the EW/TeV scale, two Higgs states cor- 
respond to the MSSM Higgs states. The s-decoupling limit in the n/sMSSM and 
UMSSM yields two CP-even Higgs bosons with similar masses and couphngs to those 
of the MSSM with one extra decoupled Higgs. The s-decoupling limit is often achieved 
in the UMSSM since the strict azz' mixing constraint must be obeyed, and requires 
s to be at the TeV scale. In this case, the mass of the decoupled Higgs scales with 
s. However, the s-decoupling limit is not always required in the UMSSM as either 
a delicate cancellation of the mixing term in the Z — Z' mass-squared matrix that 
requires tan/5 ~ \ tt^-, or ^ suppression in Qu at large tan/3 can evade the mixing 
constraint. These fine-tuning scenarios do allow s to be lower, but does not often 
result in a dramatically reduced lightest Higgs mass. In the n/sMSSM, the lightest 
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Higgs boson decouples as it has vanishing mass and is singlet dominated while H2 
and i?3 correspond to the MSSM Higgs bosons. However, the NMSSM does not have 
this behavior. Although the s-decoupling limit provides two MSSM-like Higgs bosons, 
one becomes massless at tree-level and the other scales as a/ \/2KSiXeW esc 2/3, while the 
singlet Higgs boson mass scales with k and s. This departure from the s-decoupling 
behavior of the other models is provided by the cubic self-coupling of the singlet field 
in the superpotential. 

(iv) Weak boson couplings of the Higgs bosons are generally reduced from those of the 
SM, which translates to lower Higgs production rates. However, the production rates 
can be enhanced kinematically since the Higgs mass can be lower than the SM mass 
limit. Branching fractions may be larger than in the SM due to the suppression of 
the total width if the dominant bb decay mode is suppressed by mixing effects. The 
H ^ gg partial width can be either enhanced or reduced due to both enhancements 
of couplings to fields running in the loops and their interference effects. 

(v) The branching fraction for ^ 77 can be enhanced significantly in the n/sMSSM 
and NMSSM, providing more opportunity for Higgs discovery. Interference effects aid 
the enhancement of the overall decay width, as in H ^ gg. 

(vi) Non-SM decays can become important if allowed. The hghtest Higgs boson in the 
n/sMSSM can have a high branching fraction to light neutralino pairs if kinematically 
allowed. This decay width can be as large as a few MeV and contribute significantly 
to the total width of the Higgs. However, in the n/sMSSM, much of the allowed 
kinematic region with M^o < 30 GeV may be disfavored from the prediction of a high 
Xi relic density. The Hi A2A2 decays are not favored in the UMSSM since s must 
be large to avoid the Z — Z' mixing constraint, which in turn pushes the allowed values 
of beyond the allowed region for the decay. When allowed, this decay mode is 
dominant in the NMSSM and n/sMSSM. 

(vii) The total decay width of the lightest Higgs boson can be enhanced by many orders 
of magnitude due to the large partial widths for the non-SM modes XiXi ^^id AiAi. 
Decays to Xi pairs make a contribution to the total width that is no larger than a few 
MeV, while AiAi decays lead to significant total width enhancements. 
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APPENDIX A: THE SECLUDED ?7(1)'-EXTENDED MSSM (sMSSM) 



In the UMSSM a single standard model singlet field S plays two roles: the expectation 
value of the scalar component generates an effective /^-parameter and also generates a mass 
for the Z'. There is some tension between maintaining a small enough /lefj = hs{S) and 
generating a sufficiently large Mz' and small Z — Z' mixing. As we have seen, this is most 
easily resolved by choosing a large value for s and a small hs- This typically leads to two 
CP-even MSSM-like Higgs scalars, one heavy CP-even state that is largely singlet, and one 
MSSM-like CP-odd Higgs. Similarly, if gaugino mass unification holds the UMSSM in this 
limit involves four MSSM-like neutralinos, and two heavy neutralinos which are mixtures of 
Z' gaugino and S singlino 0. 

The secluded f/(l)'-extended model (sMSSM) jl^ separates fi^s from the Z' mass by 
introducing four standard models singlets^^, S, Si, S2, S3. All of these are charged under 
U{iy and contribute to Mz', but only 5* contributes to /leff. Moreover, in an appropriate 
(decoupling) limit there is an F and D fiat direction along which 5'i^2,3 acquire large vevs, 
so that Mz' is naturally much larger than /Xes near that limit. The most general form of 
he model involves a complicated Higgs and neutralino spectrum, which was studied in Ref. 
4^. Here, we discuss the decoupling limit in which (5*1), {S2), {S3) are naturally large (TeV 
scale), and show that in that limit the >S'i_2,3 Higgs and neutralino states (and the Z' gaugino) 
approximately decouple, and that the Higgs and neutralino spectrum of the S, Hu, and Hd 
are identical to those of the nMSSM (this was commented on for the neutralinos in Ref. [3|). 

The Higgs part of the superpotential of the general sMSSM is 

W = hsSHu ■ Hd + XSS1S2S3 + W\ (Al) 

where the U{1)' charges satisfy Qh^ + Qh^ + <5s = and -Qs = Qsi = Qsz = -\Qs3, and 
W = fiiSSi + fi2SS2- The original form of the sMSSM [3l assumed /ii,2 = because small 
non-zero values reintroduce a form of the /i-problem^^. Setting /ii_2 = here would lead 
to the limit of the nMSSM with C,f = 0. We have checked that = does not alter the 
qualitative features of the nMSSM. Nevertheless, we prefer to allow /xi_2 7^ here to show 



This is motivated by concrete string constructions, which often have muhiple singlets of this kind. See, 
e.g., 0- 

Initially vanishing fii are not generated by loop corrections until soft supersymmetry breaking is turned 
on because of additional global C/(l) symmetries of the superpotential. 
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the correspondence with the general nMSSM. The Higgs potential contributions are 

Vf = %Hu- Hd + piiSi + 11282? + hi {\S\''\Hu\'' + \S\''\Hd\') 

+ \l\S^''\S2? + \XsS2S^ + + |A,5i53 + /xa^P, (A2) 

Vd = ^{\H^?-\H,?Y 

+ ^9l'(^Qs\S? + QHJHu? + QHJH,? + J2Qs^S,?^ , (A3) 

3 

Koft = ml\HX + ml\H,? + ml\S\' + Y,ml\Si\' 

i=l 

+ {AhsSH^ ■ Hd - AaA,5i5253 + h.c.) 

+ {mls.SSi + mls,SS2 + m\sAS2 + h.c). (A4) 

The soft mass-squares rrigg^ and ml^^ are required to break two unwanted global U{\)pq 
symmetries (which are also broken by /ii,2 7^ 0). The Tn?SiS2 ^^^'^ can lead to CP violation and 
can be ignored (it is useful in considering electroweak baryogenesis l48]). For = yUj = 
the potential has F and D flat directions involving large values for the secluded sector 
fields 5'i,2,3- If iT^'sSi = as well the Si tree-level potential will be unbounded below for 
appropriate values of the Si soft parameters. Turning on small values for these parameters 
(e.g., Xs ~ 0.05) leads to a natural hierarchy between the ordinary {Hu^^S) and secluded 
sector fields, with typically {Si) ~ \mSi\/Xs- 

It is clear that (jSIl), jMl), and dSH) resemble the nMSSM for the (if«,d, S) sector fields 
provided one identifies 

CsM'^ ^ mlsASi)+^lsAS2) + Xs{S^){f,r{S2)+fi2{S,)). (A5) 

However, to establish the relation to the nMSSM one must show that the mixing of these 
states with the secluded sector fields are small and also that the f/(l)' D terms in ()A3|) do 
not significantly affect the masses. 

The decoupling described above does indeed occur in the limit that Xs{Si), [fnsSi{^i)V^^ ^ 
[/ij(5'i)]^/^, hgAs, \mu,d\, \ms\, and the gaugino masses are all small compared to the f/(l)'- 

1 /2 

breaking scale Mz' ~ gv [J2Qsi\{^i)?] ■ This can be seen by explicitly examining the 
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Higgs and neutralino mass matrices, which are given for /ij = in Ref. 13|, |43| and easily 
generahzed to include fii ^ 0. The derivation is straightforward but not very enlightening, 
so we will just state the results. 

• There are nine neutralinos, consisting of 5 with a mass matrix and composition similar 
to the nMSSM, and 4 additional states that have only small mixings with the nMSSM 
sector. The latter include two heavy states of masses ~ ±M^/, which are admixtures 
of the U{1)' gaugino Z' and one linear combination of the singlinos Si, and two light 
states (the orthogonal 5*^ states). The heavy states would be maximal mixtures of 
Z' and the singlino state in the exact decoupling limit, with significant deviations 
possible away from the limit. The masses of the light singlinos and the splitting of 
the heavy pair is comparable to the electroweak scale, with the typical scale of the 
former at Xs{Si). In some cases, the LSP is actually one of the light neutralinos from 
the (approximately) decoupled sector. The small mixing effects could then lead to 
decays of the lightest neutralino from the nMSSM sector, with significant collider and 
cosmological |78| implications. 

• The CP-even Higgs sector consists of 3 nMSSM-like states, and three additional states 
that are mainly in the secluded sector 5*1^2,3- One of the latter has a large mass of 
order Mz', generated by the D terms (which yield only one large mass eigenvalue), 
while the other two may be much lighter, typically at the electroweak scale (controlled 
by [A,K{S,)f') . 

• The CP-odd Higgs states include two nMSSM-like states, and two that are mainly from 
the secluded sector (the third secluded state is the eaten f/(l)' Goldstone boson). The 

1 /2 

secluded states include one with mass controlled by [AxXs{Si)] ' and a second that 
vanishes for ^1^^ = /ij = 0. 

APPENDIX B: HIGGS MASS-SQUARED MATRICES 

We present the appropriate superpotentials and minimization conditions that are used 



0, NMSSM Q, 



to determine the tree-level Higgs mass-squared matrices in the MSSM |5 
nMSSM and the UMSSM 127 1. Additionally, radiative corrections to the mass-squared 
matrices are obtained via an effective potential whose method is outlined below. 
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1. Tree- level Higgs potential 

Prom the tree-level Higgs potential 



(Bl) 



the tree-level mass-squared matrices for the CP-even and CP-odd Higgses are obtained by 



{Ml), 



d(l)id(j)j 



d(pid(pj 



where 02, ^3) represent (0^, a), and similarly for 931,2,3- 



(B2) 



a. MSSM 

In the MSSM, the superpotential and soft breaking terms are 



/'MSSM 
''-'soft 



[ J2 iMaXaXa + Bi^H^-Hd + AthtU'Q-H^ + h.c. 

\a=l,2,3 



+ M%\Qf + M^lUf + Mp,|DP + MflLp + M'i\E 



2 I 77.12 



+ ml\Hd\'' + ml\Hu\\ 
jiving a tree-level Higgs potential contributions 



Vf 
Vd 



2| u |2 I _2 I Tj |2 



(B3) 



(B4) 



(B5) 
(B6) 
(B7) 



We replace the soft mass terms in terms of the VEVs using the potential minimization 
conditions, 





2 1 




_ 4 _ 




_ 4 _ 



1 




2 1 




2 


4 




4 



2 2 , 



(B8) 
(B9) 
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This results in the tree-level CP-even Higgs mass-squared matrix elements 



{Ml),, = 
{Ml),, = 



{Ml) 



22 



— ^VdVu - B/i, 



(BIO) 

(Bll) 

(B12) 



and tree-level CP-odd Higgs mass-squared matrix elements: 

/ A ^0 \ T-t ^d'^u 

{M'L)ij = Bn- 



ViVj 



(B13) 



where v,^, = fd.u- One of the CP-odd Higgs bosons is a massless goldstone boson that is 
absorbed by the Z boson, leaving only one physical CP-odd Higgs with tree-level mass 



(B14) 



b. NMSSM 



Superpotential and soft breaking terms: 



1 



W = KSH^ ■ Ha + -kS^ + htU'Q ■ H^, 



_ /'NMSSM 
''-'soft 



->Csoft+ ( ^KA^S^ + h.C. ) , 



where 



soft 



( J2 lMaXaXa + AshsSHu-Hd + AthtU'Q-Hu + h.c.\ 

\o=l,2,3 / 



+ ml\Hd\'^ + ml\Hy,\'^ + ml\S^^ 



(B15) 
(B16) 



(BIT) 



Potential minimization condition: 



1 

d 2 



u 2 



G^ 



hi 



4 



^2 



7:Ks' + 



1 1 

2 -'-1,2 2 -"-1,2 2 22, 

T^s = -^^s^d- l^h^V^- K S + 



+ hsKS 



V2 



f hsAs 


h^KS 


A VuS 


V V2 


^ 2 


) Vd 






VdS 


2 ) 


Vu ' 




nAf^ 




J VdVu 








(B18) 
(B19) 
(B20) 
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c. n/sMSSM 



Superpotential and soft breaking terms: 

W = hsSHu ■ Hd + ^fMIS + htWQ ■ Hu, 

Here and ^5 are mo del- dependent, dimensionless quantities. 
Potential minimization conditions: 



1 




2 1 




2 


4 




4 * 



1 



/isA^ _^ hsipMl \ VuS 



V2 



4 



2^2 , / ^s^s _^ h,ipM l\ VdS 



V2 



2 1. 2 2 1 , 2 2,/ ^^^^ V^^S^iH ^'i^'''^ 



(B21) 
(B22) 



(B23) 
(B24) 

(B25) 



d. UMSSM 



Superpotential and soft breaking terms: 

W^hjH^-Hd + htU'Q-H^, 



/'UMSSM 
''-'soft 



= -Aoft + I -MyXvXy + h.C. ] . 



(B26) 
(B27) 



We do not include the possible terms related to exotic chiral fields to cancel anomalies. 
The exotic terms depend on details of the model and we assume the masses of the exotics 
are heavy enough to be decoupled from the EW scale phenomenology that we are interested 
in. 

Potential minimization condition: 



Vd + -, 



q2 

— - hg - QHaQHuQv 



+ 
1 

2 
+ 



V2 Vd ' 



(B28) 



_ 4 
\/2 Vu ' 



hi - QHaQHudl' 



1 



- 2 [h's + Qh.Qs9i''] s' 



(B29) 



i [hi + Qh,Qs9i''] vl - \ [hi + Qn^Qsgv'] vl - \QlgyH^ + (B30) 
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2. One-loop radiative correction to the Higgs potential 



The one-loop correction of the mass-squared matrices for the CP-even and CP-odd Hig- 
gses are obtained as 



(M\) 



5.. 



1 dV^ 



5. 



1 dV^ 



Vi 



(B31) 



The second terms describes the effect of the shift of the Higgs soft mass-squared terms 
(m^, m^, ml) at the minimum so that we can still use the tree-level values of the minimization 
condition. The scalar top Yukawa couplings are expected to be dominant over the other 
model-dependent couplings. This provides a handle on the radiative corrections that is 
mo del-indep endent . 

The Coleman- Weinberg correction to the scalar potential with the stop 1-loop is given 

by 



327r2 



.i=i 



2 

m- 



3 \ 0-4/1 "^t 
, - 2m; In 

2j * V g2 



^2 



(B32) 



where Q is the renormalization scale in the DR scheme, and mf = hf \H^\^. 
The physical stop masses are given by 



H,2 = l^^^tT ^V(tr^r)2-4det^^ 



(B33) 



The stop mass-squared matrix is 



Mr 



where 



and 



Ml^ = Ml + hl\H^\ 



Ml^ = Ml + hl\H^\\ 



[MSSM], 

ht {AlHl* - hsSH^a) [NMSSM/n/sMSSM/UMSSM] . 



(B34) 

(B35) 

(B36) 
(B37) 



and is generally valid for complex At, fi and Higgs field values. However, we assume these 
parameters are real. We keep only dominant top Yukawa couplings and neglect gauge 
couplings and Yukawa couplings of other particles. 
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APPENDIX C: NEUTRALINO MASS MATRICES 



The neutralino sector is extended by the singhno, S, and in the case of the UMSSM, the 
Z'-mo, Z'. The mass matrix for the UMSSM is given by 



My-O = 



Ml 
M2 

-giVd/2 g2Vd/2 
giVu/2 -g2Vu/'2 



5'2Vd/2 -g2Vu/'2 

-/^efF 








\ 













-l^eSVu/s -l^eSVd/s gi'QsS 

gi'QHaVd gi'QHu'Vu gvQss My 



For the NMSSM, we have 



Ml 






M2 



V 



(Cl) 



(C2) 



-giVdl'^ giVul'^ 

g2Vd/'2 -5'2^^«/2 

-giVd/'2. g2Vd/'2- -/^efr -IXeSVu/s 

giVu/2 -g2Vu/'2 -fics -i^csVd/s 

-l^eSVu/s -l^eSVd/s \/2kS J 

For the n/sMSSM, the neutrahno mass matrix is the same as the NMSSM with the hmit 
K ^ taken. Note that in the s-decouphng hmit, neutrahno states in addition to those in 

the MSSM decouple with masses M^o ^ = \ (Mv ± ^/MlT^gl^Q0^ in the UMSSM, and 
M^o = ^/2ks and in the NMSSM and n/sMSSM, respectively. 



APPENDIX D: APPROXIMATIONS OF THE TREE-LEVEL HIGGS MASS IN 
VARIOUS LIMITS 

We provide approximations to the CP-even Higgs masses and exact CP-odd Higgs masses 
for the extended models considered. We replace hs with /icfr/("^) and assume /icff is at 
the EW scale in the following derivations. The approximate mass eigenvalues of the Higgs 
mass-squared matrix can be found when the mass-squared matrix elements satisfy particular 
hierarchies. 
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1. Hierarchy types 



The entries of the matrix {A, B, C) are assumed to have the same order of magnitudes 
where the powers of e <^ 1 represent the specific hierarchy. After keeping terms of up to 
(9(e^), approximate eigenvalues may be analytically given. 

Type-1 : If the dominant terms in the matrix are of the form 



A eC^ 
eC e'^B 



(Dl) 



where A is a scalar, S is a 2 x 2 matrix, and C is a 2 x 1 column vector, then the mass-squared 
matrix can be transformed to an approximate block diagonal form using 



1 - ie^rtr -ert 



er 



1 - ie^rrt 



(D2) 



where T = C/A then. 



A + e^C^C/A 0(e3) 

0(e3) e\B-CCyA) 



'11 



\ 



Ml - 



Ml, 



Ml 



^h^h (^13) 



(D3) 



This type of hierarchy is equivalent to that found in the appendix of [3| • 
Type-2 : We find this method by an analogy of the method for Type-1. 



A tC 
eCt t^B 



(D4) 



where now A is a 2 x 2 matrix, i? is a scalar, and C is a 2 x 1 column vector, we obtain the 
unitary transform on M^: 



A + Y [CC\A-^)'^ + lA-^CC^ - A-iCCt(A-i)tA) 0{e>) 

0{e-) ^{B-C^A-^C) 

0[e^) E 



(D5) 
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with 



ert i-ie2ptp 



2 

and r = A~^C. The corresponding sub-matrices are given by 



-'^-'ii TWLWl^irWL -'^-'12 2 (M?„\-^-M?,M'L 



, ,2 _ 1 Af ( A^ii - M^, ( (Mf3 )2 +(M|3)^) ^2 , M^^jMl^M^^- M^^Mj^) 

^^h2 2 ^g^F^Mp^ ^^^22+ (M^pIlMpi^ 



^ = 

»^ _ /j^2 , (M,^3)2M|,-2M,^,Mj^3M|3+M,\(M|3)^ \ /q^n 



The entries of A are assumed to be of 0{M^) + O(e^M^), and B is assumed to be of 
0{e^M^) where M^. is a common mass scale for dominant mass terms. If this hierarchy is 
not estabhshed, for example by some cancellation which disrupts the hierarchy, e.g. | (M^2)^ ~ 
-^11-^22! ^ -^c > ^^is approximation may not work. 



2. Tree-level mass expressions 

We approximate the CP-even masses in three cases (sometimes with model-specific as- 
sumptions) : 

(i) The tan/3 = 1 limit provides an exact solution. This case is meaningful since LEP 
cannot exclude this value for the extended MSSM models where the Higgs doublets 
can be mixed with a singlet. 

(ii) Both the large Mp- and large tan p hmit, where 

My = 2pLesAs CSC 2/3 (D8) 
is as the MSSM CP-odd Higgs mass-squared with /lesAg identified as /iB. 

(iii) Large s limit (s-decoupling limit). This limit provides two MSSM- like Higgs bosons 
and one singlet dominated Higgs. 

We present only dominant terms in the following CP-even expressions 
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a. MSSM (CP-even) 



The exact eigenvalues can easily be found analytically 



M^o^ = ^M| + S//csc2/3±^/Qm|-S//csc2/3) + 2M|S/xsin2/3. (D9) 



b. MSSM (CP-odd) 



Ml = 2Bncsc2(3 (DIO) 



c. Common PQ limit case (CP-even) 



(i) In the tan (3 — 1 limit, the general solutions are given by 



S2 



^ (^;Veff(^ + 2/Xeff) ± vii,s{A, - 2ii^s)Vv^ + 16s2) . (DU) 



(ii) In the large My and large tan/3 limit (My/Mc ~ 1/e and tan/3 ~ 1/e), the hierarchy 
is of Type-1 and results in 



M^o ^ -M| ± Ys^^'^zY + - cos2/3)2. (D12) 

In this limit, the lightest CP-even Higgs mass-squared always has a negative solution which 
implies this limit is not physical in the PQ limit unless loop corrections can compensate the 

negative mass-squared. 

(iii) In the large s limit {s/M^ ~ the hierarchy is of Type-2. 

+2/Xeff^ (G'2^;2 - 16/x,\ - 2Al + cos 4/3 (-G^^;^ + 2A^;))) , 

M^o ^ ^Ml + A,/Xeff CSC 2/3 ± Qm| - A,/iefr esc 2/3^ + 2M|A,/Xeff sin 2/3. (D13) 

In this limit, the hghtest Higgs mass scales as Mh^ oc which imphes that the singlet- 
dominated Higgs decouples, with mass near zero in the s-decoupling limit. Note that when 
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tan (3 ~ 1, the (MfJ^ _ Mf^M^^ term in the denominator of Eq. dHU) is of 0{e'^M^) and 
the above approximation is not vahd. However, one can take the large s hmit from the 
tan/5 = 1 exact solution of Eq. ()D11|) and arrive at the desired result. As expected, the 
other two Higgs bosons have MSSM-like masses. 



d. Common PQ limit case ( CP-odd) 



The PQ symmetry protects one of the CP-odd Higgs masses 



A, 



MX = 2/Xefr A, CSC 2/3(^1 + ^ sin2 2/5 



(D14) 



e. NMSSM case (CP-even) 



(i) In the tan/5 = 1 limit: 



M|o = M| + 2/ieffAs ^ + V2K^i,sS 



M 



HO 



^ (v'HeS {As + 2/ieff) + V2s^kA^ + As^K^ 



4s2 

± 



2\/2s^'i;\/ieff (-2/iefr + A) (16/ieff + + 2/k^ (4t;Veff (2/iefr + A^) + s^A^) 



(ii) In the large My and large tan/5 limit, the hierarchy is of Type-1. 



1/2 



(D15) 



HO 



M|o + ^ (^45/? + V2A^ ± ^ [(sM|)2 + IQv' {^l^ - M^ cos' /5)' 

/ 1 
+s^fi: f 4s\^ - 4M|sK + ^2 {2s\' - M|) + -skAI j 



(D16) 



(iii) In the large s limit: because of the k and terms, the mass-squared matrix is 
neither Type-1 nor Type-2. With an additional assumption of |k| ~ e, the leading order 
terms form a block diagonal matrix. 
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. SK 
As + — ) /ieff CSC 2P. 



As + /ieff CSC 2(3^ + 2M| (^As + fi^s sin 2/3. (D17) 

If \k\ <^ Mc/s and l^^l -C s, we obtain the large s limit of the PQ case of Eq. ()D13|1 . 
changing the leading order term of the first solution. If not, the k and terms prevent the 
Higgs sector from decoupling to two Higgs bosons with MSSM masses and very heavy or light 
singlet Higgs boson. Instead, the s-decoupling limit yields a massless Higgs, and a divergent 
Higgs {Mh2 ~ V V^KSjJicE CSC 2(3) both with ^mssm = 1 while the singlet dominated Higgs 
has a mass that scales linearly with ns. 

f. NMSSM case (CP-odd) 

The CP-odd Higgs masses are exactly given at tree-level as 
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8s2 sin 2/5 



± 



/icff^s(^^^ + 8s^) + 2V2sK/icfr(f^ + 2s^) - f VeS cos4/?(As + 2v^sk) - 6V2A«;S^k sin 2/5 

1/2 



+96s^K/ieff sin 2/5 ( 2A«;S^(a/2As + sk) - 3 V2v4^wVeff sin 2/5 



(D18) 



5. n/sMSSM case (CP-even) 



(i) In the tan/5 = 1 hmit: 



4s2 
± 



'fxls (As - 2/Xefr)' + 165^) - iV2sv^fieS {As - 2/ieff) e5M,3 ^ gs^elM,^ 



1/2 



(D19) 



(ii) In the large My and large tan/5 limit, the hierarchy is of Type-1. 



V2s 



± 



2s 



(sM| - V2^sMlf + I6v\^ls - cos^ PY + aV2MIs^sMI .(D20) 



1/2 



(iii) In the large s limit and small ^5 limit {^sM^/M^ ~ e), the hierarchy is of Type-2. 



M 



HO 



1 /ieflf sec^ 2/5 
^ 2G2A, 



(32/i^g sin 2(3 Al - G^f ^ sin=^ 2/5 [^^l^ + A^) 



+2/Xeff^ (G^^^^ - 16/i2jj - 2v42 + cos 4/5 (-G^^ + 2^2))) - 



Ml, ^ ^M| + ^/ieffCsc2/5±y QM|-^/iefrcsc2/5^ + 2M|^/Xeff sin 2/5. (D21) 

The additional assumption of a small ^s^t term provides an additional negative term in 
M^. This mass is similar to that of the PQ limit in Eq. (fm3|) when isMljMl < M^js. 
However, as s is increased, the — term is dominant and the mass diverges from the PQ 

case. This limits ^5 < to avoid negative Higgs mass-squared solutions that are unphysical. 
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However, effects from radiative corrections may alleviate this condition and allow positive, 
but small values of ^5. Along with the NMSSM, when tan/5 ~ 1, the above approximation 
is not valid, but the large s limit from the tan /? = 1 exact solution of Eq. ()D11|) may be 
taken to arrive at the desired result. 

h. n/sMSSM case (CP-odd) 



The CP-odd Higgs masses is exactly given at tree-level as 



8s^sm2(3 



fiesA.iv^ + 8s^) - AsV^fXes COS 4/3 + Sv^s/ieff^FM^^ - 4^2.3 sm2(3^sM^ 



± 



l^esAsiv^ + 8s^) - A,t;VeffCos4/3 + Sv^s/iefrCFM^ - aV2s sm2/3^sM^ 

1/2 



+64s/iefr sin 2/3 (2s{V2AsS + 2^fMI)^sM^ - V2AsV^fieS sin 2(3^fMI 



(D22) 



i. UMSSM case (CP-even) 



(i) In the tan/3 = 1 and Qh^ = Qhu limit: 



M^O = M| + 2/ieff^. - 



1 / 1 

Mlo = ^ UVeff {A, + 2/iefT) + (4s2 + (^|,g| ± [v^^lf, {As - 2/iefr)' (t;' + 16s^) 

(D23) 



1/1/ . o 2X^4 ,o4 



(20s2 - /ieff (A, - 2/iefr) (^I^QI + (4s' + t;') g%Qs 
(ii) In the large My and large tan/3 limit, the hierarchy is of Type-1. 



1/2^ 



1 



1 



fro ^ -M7 + -Ml, ± 

^ 2^2^ 2s 



s'Ml + s'Ml, + 16v' (^,\ - cos' /3)^ + s'g',, {v'QI^ + s^Qlf 

(D24) 



+-s\^gl, {GWQI.^ - Gh^Ql + 32 - cos' /3) Qh^Qs) 

(iii) In the large s limit: if we change the basis of the matrix from {H^-, H'^, S} to 
{5*, H^, H'^}, the hierarchy is of Type-1. 



M 



Mi, 



Z'l 



Mho ~ + Asfies CSC 2/3 ± W Qm| - A./ieg esc 2/3^ + 2M|A,/ieff sin 2/3. (D25) 
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In this limit, one of the Higgs masses scales as Mh ~ Mz' which scales as s in the s- 
decoupling limit. This yields a heavy singlet dominated Higgs which decouples while the 
other two CP-even states become MSSM-like. 

j. UMSSM (CP- odd) 

The CP-odd Higgs mass in the UMSSM is equivalent to that in the common PQ limit 
(e.g. Eq. dnil) 

APPENDIX E: ADDITIONAL ONE PARAMETER PLOT RESULTS 

The one parameter plots performed in section |3 show the lightest CP-even and CP-odd 
Higgs masses vs. tan/? and s. In Fig. ^1 we show the complete Higgs mass spectrum of 
each model and the Peccei-Quinn limit versus tan jS and s. For comparison, the Higgs mass 
spectra of the mSUGRA/constrained-MSSM model can be found in Refs jSOf. 

In this figure, the level crossings between Higgs states are apparent as s and tan/3 are 
varied. The general behavior of the heaviest CP-odd mass-squares are similar as they scale 
as tan/3 and s with a minimum near tan/? ~ 1. The first and second CP-even Higgs states 
cross at intermediate tan j3 in all the models for these set of parameters, a feature not present 
in the MSSM. This suggests that the lightest CP-even Higgs typically has a maximal mass 
at some intermediate tan/3. 

The s plots shows a relatively constant lightest CP-even Higgs mass with minor variations 
when level crossings occur. Since the mass-squared difference between the charged Higgs 
and heaviest CP-odd Higgs is, excluding model-dependent variations, on the order of 
both masses have the same dependence on s and tan f3. Departures from the similarities 
in mass signify a model dependence in Mfj± and M^2- For large s, the behavior of the 
charged Higgs boson is effectively the same among the models as its mass scales with s, 
see, e.g., Eq. fH^ . The heaviest CP-even and CP-odd Higgs bosons masses also scale with 
increasing s for most of the models in this parameter range. However, the corresponding 
n/sMSSM masses decrease for increasing s at low s. This occurs when Aiss becomes the 
dominant element of the mass-squared matrix. This occurs for all the models at low enough 
s (typically beyond the lower limit of s allowed by the theoretical constraints). However, the 
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MSSM NMSSM n/sMSSM 




Ian p Ian P Ian p 



UMSSM All models, PQ limit 




FIG. 15: Higgs masses vs. tan/3 and s for the MSSM (with fi = hss/V2), NMSSM, n/sMSSM, 
UMSSM, and the PQ hmit for the extended models. Only the theoretical constraints are applied. 
Input parameters of = 500 GeV, At = I TeV, Mq = M^ = l TeV, k = 0.5, A^ = -250 GeV, 
Mn = 500 GeV, = —0.1, ^5 = —0.1, hg = 0.5, 9e& = — tan~^ ^J^, and the renormalization scale 
Q = 300 GeV are used. Note that the U{1)pq symmetry allows only one CP-odd Higgs boson to 
be massive. 



term specific to the n/sMSSM shifts the point where this cross-over from Mfj^ ~ 
to Mjj^ ~ ^ occurs to higher values of s, resulting in the distinctive departure from the 
behavior of the other models at low s after theoretical constraints are applied. 
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